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4  INTRODUCTION: 

Aberrant  expression  of  receptor  or  cytosolic  tyrosine  kinase  genes  and,  in  particular,  their  hyper¬ 
expression  are  common  phenomena  in  breast  cancer,  which  are  believed  to  alter  cell  growth  and 
response  to  external  signals  such  as  growth  factors,  hormones  etc.  Knowledge  about  the  relative 
levels  of  expression  of  many  tyrosine  kinase  genes,  all  at  the  same  time,  might  contribute 
significantly  to  a  better  understanding  of  the  processes  of  tumor  development  and  progression.  We 
are  developing  a  rapid  assay  that  will  use  innovative  cDNA  micro-arrays  carrying  small  amounts  of 
individual  tyrosine  kinase  gene-specific  targets  to  simultaneously  determine  the  expression  levels  of 
up  to  100  tyrosine  kinase  genes  using  a  small  number  of  cells.  Three  years  of  research  and 
development  will  lead  to  discovery  of  a  set  of  gene-specific  markers  associated  with  breast  cancer 
progression  and  a  simple  device  capable  of  performing  inexpensive  expression  profiling  of  these 
markers.  The  research  and  development  efforts  in  the  first  year  of  this  3-year  project  focused  on  the 
esign  and  testing  of  robotic  instruments  to  prepare  DNA  micro-arrays  and  the  preparation  of 
prototype  arrays  carrying  sets  of  more  than  50  gene-specific  tyrosine  kinase  fragments.  The  second 
year  effort  (i.e.,  the  reporting  period  1  Aug  00  -31  Jul  01)  was  directed  towards  the  optimization  of 
cDNA  preparation,  labeling,  hybridization  and  detection  protocols  as  well  as  the  molecular  cloning 
and  sequencing  of  breast  cancer- specific  tyrosine  kinase  gene  transcripts. 

BODY: 

Here,  we  report  our  progress  as  it  relates  to  the  approved  ’Statement  of  Work’. 

Task  1.  Identify  tyrosine  kinase  (tk)  genes  expressed  in  normal  and  neoplastic  breast  tissues 
Prepare  cDNAs  front  ten  cell  lines  and  ten  frozen  tissue  specimens  (months  1-18) 

We  isolated  total  RNA  from  16  different  breast  cell  lines  and  ten  frozen  tissue  specimens 

Exponentially  growing  cell  lines  or  isolated  RNA’s  were  provided  to  us  by  the  following 
collaborators: 

Paul  Yaswen,  Ph.D„  Martha  Stampfer,  Ph.D.,  Daniel  Callahan,  Ph.D.,  and  Ruth  Lupo,  Ph.D.,  from 
LBNL,  and  Chris  Benz,  M.D.,  from  the  Buck  Institute  for  Aging  Research,  Novato,  CA.  The  actual 
number  of  cell  lines  or  RNA’s  provided  is  higher,  because  we  processed  several  aliquots  of  what  is 
believed  the  same  cell  line.  For  example,  we  received  MCF-7  cells  from  D.  Callahan  and  MCF-7 
RNA  from  P.  Yaswen  and  C.  Benz. 


Table  I.  Breast  cancer  cell  lines  used  in  our  experiments. 

ZR75  - 

BT474 
BT549 
184A1 
184B5 
T47D 


MDA-MB-436  plus  Vit.  A 
600MPE 


SKBR3 

BT468 

MCF-7 

184A1TH-6 

1 84B5-ME 

MDA-MB-453 

MDA-MB-436  plus  pCMr3.1 

HMEC 
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Normal  breast  tissue  and  breast  cancer  tissue  specimens  were  obtained  from  the  University  of 
California  Comprehensive  Cancer  Center  (Joe  W.  Gray,  Ph.D.,  B.M.  Ljung,  M.D.  and  K.  Chew). 
We  received  a  total  of  10  frozen  tissue  specimens  representing  3  normal  tissues  and  7  cancer  tissues. 
RNA  was  extracted  using  a  commercial  kit  (Qiagen)  and  transcribed  into  cDNA  immediately  after 
isolation.  Remaining  RNA  was  stored  at  —80  degrees. 


We  prepared  cDNAs  from  the  RNA  by  random  priming  and  reverse  transcription.  The  names  of 
e  various  cell  lines  are  shown  in  Table  I.  Commercial  kits  (Qiagen,  Roche,  Ambion)  were  used 

for  all  steps.  Typically,  1  pg  of  total  RNA  produced  sufficient  quantities  of  cDNA  for  cloning 
and/ or  repeated  micro-array  analyses. 


1.2  Perform  RT-PCR  reactions  and  clone  PCR  products  in  plasmids  (months  3-18) 

We  used  about  100  ng  of  DNA  in  PCR  reactions  to  amplify  tk-specific  cDNA  fragments  with 
our  mixed-base  F-/R-TYRK  primers.  PCR  amplification  was  performed  for  35  cycles  using  mixed 
base  primers  that  bind  to  the  conserved  sequences  of  domains  VII  and  IX  of  the  tyrosine  kinase 
genes  (F-TYRK:  5  '-GGGCGTCAGAARRTNRSNGAYTTYGG-3 R-TYRK:  5'-GCGCGGGCC- 
CRWANSHCCANACRTCNSA-3 ').  Each  cycle  consisted  of  a  denaturation  step  of  30  sec  at  94°C, 
primer  annealing  at  53°C  for  60  sec  and  primer  extension  for  120  sec  at  72°C.  Most  PCR  products 
appeared  as  a  single  broad  band  of  the  expected  size  (-160-170  bp).  The  amplification  products 
were  separated  on  a  4%  agarose  gel,  and  a  slice  containing  fragments  of  approximately  160-180  bp 
was  excised.  The  gel  slice  was  rinsed  with  water  and  melted  by  heating  in  1 00  pi  water.  About  2  pi 
of  this  solution  was  transferred  to  200  pi  of  PCR  buffer  containing  the  modified  F-TYRK/R-TYRK 
primers  designed  to  include  deoxy-UMP  residues  and  a  Not  I  restriction  site  at  their  respective  5'- 
ends  (F-TYRKU:  5  '-CUACUACUACUAGCGGCCGCAARRTNRSNGA YTTYGG-3  '•  R-TYRK- 
5  '-CAUCAUCAUCAUGCGGCCGCCCRWANSHCCANACRTCNS-3 ').  This  gave  PCR  pro¬ 
ducts  of  the  expected  size  (about  190  bp)  suitable  for  treatment  with  Uracil  DNA  Glycosylase 
(UDG)  and  cloning  into  pAMPl  (Gibco/LTI).  Transformation  was  performed  using  MAX 
fficiency  DH5a  competent  cells  (Gibco/LTI).  Cells  were  incubated  overnight  at  37°C  on  LB  plates 
containing  lOOpg/ml  ampicillin,  50pg/ml  X-gal,  and  ImM  IPTG  (Gibco/LTI).  Ampicillin  resistant 
clones  were  picked  from  agar  plates  and  their  insert  sizes  were  determined  by  agarose  gel  analysis 

of  PCR  products  generated  with  vector-specific  PCR  primers. 

1.3  Perform  pre-screening  with  known  tk fragments,  cDNA  sequencing  and  database  searches 
(months  5-22) 

The  DNA  from  plasmid  clones  with  inserts  of  about  125-190  bp  was  isolated,  fingerprinted  or 
screened  against  known  tk  genes,  and  ‘novel’  clones  were  sequenced  at  the  UC  Berkeley, 
Biochemical  Core  Facility.  The  list  of  genes  cloned  and  identified  so  far  contained  several 
sequences  reported  to  have  transforming  activity,  such  as  trk,  axl/ufo  or  to  be  overexpressed  in 
various  types  of  cancer.  Our  present  panel  of  kinase  gene  tags  used  to  prepare  cDNA  micro-arrays 
contains  60  genes  among  them  four  novel  sequences  and  HLA-A,  which  happens  to  be  amplified  by 
our  PCR  primers.  In  the  reporting  period,  we  noticed  an  elevated  fraction  of  clones  that  were 
sequenced  and  revealed  a  previously  sequenced  insert,  suggesting  that  our  pre-screening  procedure 
was  not  sufficiently  stringent.  For  example,  one  batch  sequencing  procedure  of  20  plasmid  clones 
showed  previously  identified  tk  fragments  in  18  out  of  20  reaction,  while  the  remaining  2  clones 
contained  novel  sequences  without  matches  in  the  Genbank  database. 
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1.4  Add  novel  clones  to  the  panel  of  expressed  tk  gene  fragments  (months  6-24) 

As  of  June  2001,  we  cloned  and  partially  characterized  about  240  tk  fragment  containing  plasmid 
clones  derived  from  breast  cancer  tissues  and  cell  lines.  These  clones  join  more  than  500  tk 
fragments  containing  clones  that  were  previously  isolated  from  thyroid  and  prostate  tumors.  From 
these  clones,  DNA  was  isolated,  bound  to  nylon  filters  and  prescreened  with  probes  prepared  from 
known  tk  clones.  Following  this  pre-screening  step,  we  performed  cDNA  sequencing  and  database 
searches.  In  about  100  tk  containing  plasmids  that  were  sequenced,  we  found  two  potentially  novel 
tk  genes.  The  other  sequenced  clones  contained  two  known  kinases  that  we  added  to  our  panel.  The 
remaining  clones  contained  sequences  that  were  already  part  of  the  panel  suggesting  insufficioent 
stringency  during  the  pre-screening  process.  At  this  point  we  decided  to  focus  our  efforts  on  the 
optimization  of  the  PCR  amplification  (see  2.3,  below).  And  added  novel  clones  to  the  panel  of 
expressed  tk  gene  fragments.  DNA  isolated  from  this  panel  of  tk  gene  fragments  is  then  arrayed  and 

printed  onto  glass  slides  to  measure  tk  gene  expression  in  cell  lines,  normal  tissues,  and  tumor 
tissues. 

We  initiated  a  search  for  full  length  cDNA  clones  for  some  of  the  novel  tk  sequences  found 
expressed  in  breast  cancer.  We  will  hire  a  local  company  (Pangene  of  Fremont,  CA)  to  screen  the 
proprietary  cDNA  libraries  for  full  length  clones. 

Task  2.  Measure  tk  gene  expression  in  cell  lines,  normal  and  tumor  tissues 

2.1  Prepare  DNA  microarrays  carrying  about  100  different  sequences  (months  3-24) 

"^^fabrication.  tyrosine  kinase  (tk)  cDNA  arrays  on  glass  slides  was  done  in  house  using 
the  laboratory-built  robotic  system  to  print  cDNA  micro-  arrays  with  about  124  target  spots  per 
array.  The  arrayer  can  print  up  to  91  slides  in  a  single  run.  Routinely,  we  print  approximately  50 
slides  per  run,  with  duplicate  arrays  on  each  slide.  The  stainless  steel  printing  pin  has  an  open  slit 
in  the  tip.  Each  loading  takes  up  approximately  1  pi  of  DNA  solution  and  allows  continuous 
printing  of  more  than  100  dots.  A  full  384-well  plate  of  DNA  samples  can  be  arrayed  onto  50 
slides  (in  duplicate)  in  8  hrs,  including  the  time  needed  for  repeated  wash  and  drying  steps  in 
each  printing  cycle.  We  concentrated  our  PCR  products  and  dissolved  DNA  in  50%  dimethyl 
sulfoxide  (DMSO)  at  a  concentration  of  250ng/pl.  Spot  diameters  on  poly-L-lysine  coated  slides 
ranged  from  80-125  pm.  The  arrayer  has  a  positioning  precision  of  approximately  15pm, 
therefore,  an  array  with  200pm  center-to-center  distance  between  spots  is  sufficiently  spaced. 

Standard  microscope  glass  slides  were  cleaned  thoroughly  with  concentrated  sodium  hydroxide 
and  ethanol  before  being  coated  with  poly-L-lysine.  These  steps  provided  a  positively  charged  layer 
to  bind  DNA  to  the  slides.  The  coated  slides  were  stored  in  ambient  air  in  closed  slide  boxes  two 
weeks  before  they  were  used  for  printing.  The  locations  of  arrays  were  marked  on  the  reverse  side 
of  the  slides  with  diamond  marker  pen.  The  slides  were  placed  upside  down  in  a  chamber  filled  with 
l  x  SSC  and  re-hydrated  for  5  to  1 5  minutes  (depending  on  the  size  of  the  array).  Spots  are  not 
allowed  to  swell  too  much  and  run  into  each  other.  Hydrated  slides  were  dried  on  a  hot  plate  set  to 
70  C  for  3  seconds.  These  slides  were  then  UV  cross-linked  using  a  Stratalinker  1800  (Stratagene) 
set  to  65  mJ.  Following  cross-linking,  slides  were  incubated  in  blocking  solution  (0.2 1M  succinic 
anhydride  m  l-methyl-2-pyrrolidinone,  0.06M  sodium  borate  for  20  minutes  (pH  8.0)).  The  slides 
were  then  denatured  in  boiling  water  for  2  min,  and  submerged  in  95%  ethanol  for  1  minute.  After  a 
centrifugation  at  500  rpm  for  5  min  to  remove  traces  of  ethanol  or  water,  the  arrays  were  ready  for 
hybridization  or  storage. 

Quality  control  for  the  tk  arrays  was  monitored  throughout  the  tk  array  fabrication  process.  All 
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reagents  used  for  printing  and  coating  of  slides  were  prepared  fresh  and  filtered.  Before  post¬ 
printing  processing  of  the  arrays,  morphology  of  the  spot  and  DNA  distribution  within  the  spots 
were  examined  with  DAPI  staining  and  visualized  in  a  Zeiss  fluorescence  microscope.  After  post- 
pnnting  processing,  arrays  were  stained  with  POPO-3  dye  (0.01  mM)  and  scanned  at  Cy3  channel 
to  ensure  the  retention  of  DNA  on  the  arrays  for  hybridization. 


2.2  Optimize  hybridization  conditions  to  provide  quantitative  information  (months  6-24) 


We  isolated  total  RNA  from  breast  cancer  cell  lines  or  frozen  tissues  sections.  We  then  prepared 
cDNAs  from  the  RNA  by  random  priming  and  reverse  transcription  using  commercial  kits  (Qiagen, 
Roche,  Ambion).  Typically  1  pg  of  total  RNA  can  produce  sufficient  quantities  of  probe  for  all 
microarray  analyses.  We  used  about  100  ng  of  cDNA  in  PCR  reactions  to  amplify  tk-specific 
cDNA  fragments  with  our  mixed-base  F-/R-TYRK  primers.  Following  30  cycles  of  PCR,  the 
amplification  products  were  separated  on  a  4%  agarose  gel  stained  with  ethidium  bromide,  and  a 
slice  that  contained  fragments  of  approximately  160-180  bp  was  excised.  The  gel  slice  was  rinsed 
with  water  and  melted  by  heating  in  1 00  pi  water.  About  2  pi  of  this  solution  was  transferred  to  200 
pi  of  PCR  buffer  containing  the  F-/R-TYRK  primers.  The  tk  fragments  were  then  amplified  a 
second  time  and  the  products  concentrated  via  isopropanol  precipitation.  The  tk  fragments  (about 
400-600  ng  per  reaction  (Hsieh  et  al.  2001))  were  then  labeled  with  either  Cy3-dUTP  (red 
fluorescence,  shown  green  in  the  figures)  or  Cy5-dUTP  (infrared  fluorescence,  shown  red  in  the 
figures)  by  random  priming  (Life  Technologies). 

We  optimized  hybridization  and  wash  conditions  to  obtain  specificity  and  reduce  background 
fluorescence  to  the  level  in  which  quantitative  information  could  be  obtained  reproducibly.  Slides 
were  placed  inside  hybridization  chambers  with  a  cover  fastened  tightly  by  clamps.  Multiple  drops 
of  3X  SSC  were  added  around  the  arrays  on  the  slide  to  maintain  the  humidity  and  prevent  diyingof 
hybridization  solution.  Hybridization  was  then  performed  under  glass  coverslips  overnight  at  65  °C. 
Next  day,  the  arrays  were  then  washed  in  3  changes  of  wash  solution  (lx  SSC+  0.03%  SDS)  at 

ln°o^tfmperatUre  for  5  minutes’  °-2x  ssc  at  room  temperature  for  1  minute,  and  then  0.2x  SSC  at 
50  C  for  1  minute.  After  the  washes,  the  slides  were  dried  immediately  by  centrifugation  at  500 
rPm]  *or  ^  .n^n  before  scanning.  Current  sample  preparation  and  hybridization  procedures  have 
produced  signal  to  background  (noise)  ratios  of  more  than  100  (in  the  green  channel)  to  more  than 
200  (m  the  red  channel). 

The  performance  of  the  system  was  tested  by  hybridization  of  fluorochrome-labeled  tk  gene- 
specific  PCR  fragments  onto  our  tk-specific  DNA  micro-arrays.  In  preliminaiy  experiments,  we 
used  a  first  generation  of  DNA  micro-arrays  comprised  of  48  kinase  genes  previously  identified  to 
be  expressed  in  thyroid  tumors  (not  shown).  This  allowed  us  to  optimize  hybridization  and  wash 
conditions.  We  then  prepared  DNA  micro-arrays  carrying  a  more  extensive  panel  of  tk  genes  plus  a 
few  control  genes  (Table  n.  We  found  that  hybridization  to  these  arrays  (Figure  1)  provided 
reproducible  information  relevant  to  tumor  progression  and  differences  between  cell  lines.  For 
some  of  our  experiments,  we  want  to  compare  the  expression  levels  of  tk  genes  in  a  tissue  sample 
against  a  constant  standard  rather  than  against  a  second  tissue  sample.  Here,  the  small  size  of  our 
array  is  an  advantage.  We  can  readily  create  a  reference  standard  to  compare  tk  expression  in  a 
given  tissue  simply  by  making  a  cocktail  containing  equimolar  amounts  of  each  tk  spot  on  the  array, 
uus  artificial  cocktail  mix  of  DNAs  is  then  labeled  by  random  priming  with  Cy3-dUTP  or  Cy3-dUTP 
in  an  identical  manner  as  the  cDNAs  derived  from  mRNA. 
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2.3  Optimize  PCR  parameters  for  quantitative  amplification  of  target  genes  (months  12-18) 

We  have  worked  out  an  effective  method  to  generate  probe  for  our  microarray  experiments.  In 
our  original  method,  we  isolate  total  RNA  from  frozen  tissues  cell  lines  and  prepared  cDNAs  from 
the  RNA  by  random  priming  and  reverse  transcription,  before  we  amplify  tk-specific  cDNA 
fragments  with  our  mixed-base  F-/R-TYRK  primers. 

In  the  course  of  our  experiments,  we  have  discovered  that  this  probe  synthesis  method  can 
generate  quality  microarray  data.  This  method  is  particularly  good  at  assaying  the  expression  of  tk 
genes  that  are  poorly  or  moderately  expressed.  We  subjected  our  results  to  independent  validations 
with  colleagues  expert  in  the  expression  of  tk  genes  in  the  breast  cells.  We  have  the  learned  that  our 
original  method  seems  to  recapitulate  the  results  of  other  investigators  (using  Northern  blots)  for 
moderate  and  low  expressing  tk  genes.  However,  we  have  also  learned  that  highly  expressed  tk 
gene  expression  is  less  accurately  represented  in  our  experiments  than  we  would  prefer.  The  rank 
ordering  of  the  high  expressing  genes  is  maintained  (i.e.  for  tk  genes  that  are  highly  expressed,  a 
lower-expressed  tk  gene  indeed  appears  to  be  expressed  less  than  a  higher-expressed  tk  genes). 
However,  the  apparent  difference  between  highly  expressed  genes  are  less  than  one  might  expect. 

In  essence,  our  original  probe  synthesis  scheme  is  not  completely  representative  of  total  tk 
expression. 

We  believe  that  the  cause  of  this  artifact  is  due  to  the  depletion  of  degenerate  primers  in  the  early 
steps  of  PCR  synthesis.  Our  degenerate  primer  mixes  contain  equimolar  concentrations  of  each  of 
die  hundreds  of  different  primers.  Of  course,  the  population  of  tk  RNAs  within  and  tissues  sample 
is  likely  to  be  anything  but  equimolar.  It  is  probable  that  the  most  highly  expressing  genes  are 
rapidly  depleting  the  population  of  degenerate  primers  that  would  be  perfect  matches  for  that 
species.  Once  perfect  primers  are  depleted,  sub-optimal  primers  then  bind  with  baspair  mismatches 
leading  to  degrade  PCR  performance. 

We  are  currently  testing  revised  methods  to  eliminate  this  problem.  Our  new  PCR  scheme 
utilizes  new  degenerate  tk  primers  that  have  an  unrelated  18-25  bp  adapter  sequence  at  its  5’-end 
and  the  tk  -specific  sequence  on  its  3 ’-end.  The  previous  work  of  Eberwine  (Eberwine,  1996)  with  a 
T7  adaptor  sequence  for  in  vitro  transcription  suggested  that  Ml 3  or  T7/T3  primer  sequences  on  the 
5 ’-ends  should  work  well.  The  new  degenerate  primers  are  used  for  only  a  handful  of  PCR  cycles 
with  the  cDNAs  to  initiate  the  amplification  of  the  tk  sequences.  Then  the  degenerate  primers  are 
removed  with  a  spin  column,  and  the  adapter  primers  (representing  the  non-degenerate  5’-ends)  are 
added  to  the  products  of  the  first  PCR  reaction.  The  adapter  primers  then  complete  the  reaction.  By 
allowing  the  degenerate  primers  access  to  our  cDNA  for  only  a  very  few  cycles,  we  believe  we  can 
make  our  microarray  results  more  representative.  This  will  be  even  more  important  as  we  move  to 
smaller  samples  sizes.  We  have  already  tested  different  adaptor  sequences,  and  are  now  in  the 

process  to  compare  the  complexity  of  PCR  products  generated  with  either  to  old  or  the  new 
amplification  strategy. 
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Figure  1:  Results  of  hybridization  of  CY3 -labeled  tk  DNA  fragments  from  breast  cancer  cell  line 
184A1  (shown  in  green)  in  combination  with  CY5  -labeled  TK  fragments  from  normal  breast 
tissue  (shown  in  red)  reveals  specific  differences  in  the  expression  of  tk  genes. 
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Table  II:  The  organization  of  our  new  tk  DNA  microarrays  used  in  studies  of  gene  expression  in 
breast  cancer  cell  lines  and  tissue  specimens. 
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2*4  Develop  algorithms  for  array  readout  and  comparisons  between  measurements  (months  9-24) 

A  Axon  GenePix  4000  (Axon  Inc.)  array  scanner  used  to  acquire  all  images  has  a  preview 
resolution  of  40  (im  and  a  scanning  resolution  of  10  pm.  The  photomultiplier  sensitivity  can  be 
adjusted  by  the  user  during  preview  to  optimize  the  signal  intensity.  GenePix  3.0  (Axon  Inc.),  an 
image  acquisition  and  analysis  software,  was  used  to  analyze  the  images  acquired  from  our  tk  arrays 
and  provided  numerical  data  that  was  imported  into  spreadsheets  (Microsoft  Excel)  for  further 
analysis.  For  display  purposes,  the  images  were  saved  in  standard  formats  and  imported  into 
common  graphics  programs  such  as  Adobe  Photoshop.  Given  the  relatively  small  numbers  of 
arrays  that  have  been  hybridized  so  far,  it  was  sufficient  to  normalize  the  data  and  compare  numeric 
values  in  spreadsheets. 

Task  3 .  Validate  assays  for  multigene  expression  profiling  in  small  amounts  of  tissue 

3*1  Develop  software  for  databasing ,  automated  analysis  of  expression  profile  datasets  and  their 
annotation  (months  18-24) 

The  GenePix  software  provides  the  results  as  numeric  data  that  can  be  inported  into 
spreadsheets  or  databases.  We  installed  Microsoft  Access  software  which  will  be  sufficient  to 
handle  most  or  all  of  our  databasing  needs.  We  also  installed  BRB  ArrayTools  Version  1 .03 
recently  released  by  the  Biometric  Research  Branch  of  the  Division  of  Cancer  Treatment  and 
Diagnosis  at  NCI.  This  package  developed  by  Richard  Simon  and  Amy  Peng  allows  for  a  more 
comprehensive  analysis  of  DNA  microarrays  and  their  annotation.  The  software  is  presently  being 
tested  using  the  datasets  obtained  with  our  second  generation  tk  cDNA  arrays. 
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KEY  RESEARCH  ACCOMPLISHMENTS: 


•  Finished  the  isolation  of  RNA  and  preparation  of  cDNA  from  16  breast  cancer  cell  lines  and 
10  frozen  tissue  specimens 

•  Completed  the  PCR-amplification  of  tk-specific  DNA  fragments  and  cloned  the  products  into 
plasmids 

•  Pre-screening  more  than  400  breast  cancer  cell  line-derived  clones  and  sequenced  an  additional  100 
clones,  database  searches  identified  two  clones  containing  potentially  novel  tyrosine  kinase  genes 

•  Expanded  the  panel  of  tyrosine  kinase  genes  used  for  expression  profiling  and  printed  second 
generation  cDNA  micro-arrays 

•  Generated  artificial  mixtures  of  tk  DNA  fragments  to  be  used  as  reference  DNA 

•  Reconfirmed  tk  gene  expression  changes  as  breast  epithelial  cells  become  tumorigenic  and  grow 
anchorage-independent  using  second  generation  tk  micro-arrays 

•  Initiated  search  for  full  length  cDNA  clones  for  novel  tk  genes 
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CONCLUSIONS: 

This  3-year  IDEA  project  is  well  on  track  and,  has  met  its  initial  milestones.  The  soft-  and  hardware 
components  necessary  for  these  studies  were  put  in  place  in  the  first  year.  The  results  obtained  with 
RNA  isolated  from  cell  lines  and  breast  tissues  have  proven  the  hypothesis  that  changes  in  tk  gene 
expression  can  be  monitored  by  a  combination  of  PCR  using  tk  gene  family-specific  primers  and  DNA 
micro-arrays.  While  the  hybridization  to  the  DNA  micro-array  appears  to  possess  the  required 
specificity,  second  year  research  addressed  the  issues  of  hybridization  background  reduction  and 
definition  of  a  suitable  reference  DNA  probe.  Comparison  of  the  cDNA  micro-array  data  with  those 
obtained  by  Southern  blot  analyses  suggests  a  non-homogeneous  amplification  of  tk  fragments.  This  is 
now  addressed  by  an  altered  PCR  protocol  involving  new  primers.  Concordant  with  the  timeline 
presented  in  the  original  proposal,  research  and  development  in  the  third  year  will  focus  on  application 
of  the  technology  to  small  numbers  of  cells. 
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Monitoring  Signal  Transduction  in  Cancer:  cDNA  Microarray  for 
Semiquantitative  Analysis 

H.-Ben  Hsieh,  Robert  A.  Lersch,  Daniel  E.  Callahan,  Simon  Hayward,  Mariwil  Wong, 

Orlo  H.  Clark,  and  Heinz-Ulrich  G.  Weier 
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Departments  of  Urology  (SH)  and  Surgery  (MW,OHC),  University  of  California-San  Francisco,  San  Francisco,  California 

SUMMARY  This  study  targeted  the  development  of  a  novel  microarray  tool  to  allow 
rapid  determination  of  the  expression  levels  of  58  different  tyrosine  kinase  (tk)  genes  in 
small  tumor  samples.  The  goals  were  to  define  a  reference  probe  for  multi-sample  compar¬ 
ison  and  to  investigate  the  variability  and  reproducibility  of  the  image  acquisition  and  RT- 
PCR  procedures.  The  small  number  of  tk  genes  on  our  arrays  enabled  us  to  define  a  refer¬ 
ence  probe  by  artificially  mixing  all  genes  on  the  arrays.  Such  a  probe  provided  contrast 
reference  for  comparative  hybridization  of  control  and  sample  DNA  and  enabled  cross¬ 
comparison  of  more  than  two  samples  against  one  another.  Comparison  of  signals  generated 
from  multiple  scanning  eliminated  the  concern  of  photo  bleaching  and  scanner  intrinsic 
noise.  Tests  performed  with  breast,  thyroid,  and  prostate  cancer  samples  yielded  distinctive 
patterns  and  suggest  the  feasibility  of  our  approach.  Repeated  experiments  indicated  re¬ 
producibility  of  such  arrays.  Up-  or  downregulated  genes  identified  by  this  rapid  screening 
are  now  being  investigated  with  techniques  such  as  in  situ  hybridization. 

(J  Histochem  Cytochem  49:1057-1058,  2001) 
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cDNA  microarray 
expression  profiling 
tyrosine  kinase 
reference  probe 


Our  DNA  microarray  experiments  followed  the 
Brown  laboratory’s  protocols  (http://cmgm. Stanford, 
edu/pbrown/mguide/index.html)  with  modifications 
on  microarrayer  construction  and  the  following:  spot¬ 
ted  DNAs  are  homogeneous  (270-290  bp);  RNAs  were 
amplified  by  RT-PCR  with  mixed-base  tk  primers,  size- 
selected  by  electrophoresis  (180-210  bp),  and  labeled 
by  random  prime  reaction.  Typical  cDNA  microarray 
experiments  utilize  two-color  hybridization  in  which 
Cy5-labeled  sample  is  hybridized  against  Cy3-labeled 
control  to  the  same  arrays  (Schena  1996).  This  provides 
a  convenient  means  to  look  at  two  treatment  effects, 
such  as  two  time  points  (Iyer  et  al.  1999)  or  tumor  vs 
normal  tissues  (DeRisi  et  al.  1996).  However,  when 
multiple  samples  are  to  be  compared  against  one  an¬ 
other,  pair-wise  comparisons  are  inefficient.  To  address 
such  needs  and  reduce  the  complexity,  a  “standard”  or 
“reference”  probe  is  preferred.  Because  the  hybridiza- 
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tion  is  competitive  with  limited  labeled  DNA  (Cheung 
et  al.  1999),  a  pooled  mixture  that  represents  all  genes 
at  low  (but  detectable)  levels  will  be  an  ideal  reference. 
For  large-format  cDNA  microarray,  pooled  mixtures  of 
several  normal  tissues  or  cell  lines  were  used  (Celis  et  al. 
2000).  Such  pools  of  RNA  can  vary  from  preparation 
to  preparation,  and  highly  expressed  genes  can  skew  or 
saturate  the  arrays.  With  a  small  number  (58)  of  genes 
on  our  tk  arrays,  we  were  able  to  artificially  mix  indi¬ 
vidual  genes  proportionally  and  to  amplify  them  in 
large  quantity  for  a  particular  project.  Pair-wise  hybrid¬ 
izations  with  this  reference  probe  have  yielded  similar 
results  that  proved  this  approach  feasible. 

To  test  the  extent  of  the  scanner’s  intrinsic  noise, 
slides  were  scanned  12^ times  repeatedly  at  full  laser 
power.  Four  separate  pixels  belonging  to  four  inten¬ 
sity  groups  were  selected  and  their  mean,  median,  SD, 
and  CV  were  determined  (Table  1).  Results  indicated 
an  average  of  5.0%  and  6.0%  CV  for  Cy5  and  Cy3 
channels,  respectively.  Noises  were  negligible  and  did 
not  correlate  with  PMT  voltages;  the  deviations  of 
measured  intensities  were  within  experimental  errors 
(<<  10% ).  Repeated  scanning  up  to  24  times  did  not 
decrease  the  fluorescence  intensity,  eliminating  the 
concern  of  photo  bleaching. 
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Table  1  Repeated  scanning  of  the  same  array  12  times  at  full  laser  power.  Results  for  four  different  pixels  in  Cy5  and  Cy3  channels 
indicated  negligible  photo  bleaching  and  «10%  coefficients  of  variation 


Pixel  # 

1 

2 

3 

4 

Coordinate 

X  =  12390 

Y  «  9420 

Cy5 

Cy3 

X  =  12600 

Y  =  9420 

Cy5 

Cy3 

X  =  12790 

Y  =  9420 

Cy5 

Cy3 

X  —  12980 

Y  =  9630 

Cy5 

Cy3 

50%  PMT 

Mean 

1027.4 

3545.9 

416.7 

702.3 

1966.3 

2112.2 

2215.5 

1257.6 

Median 

1024.5 

3534.5 

409.5 

708.5 

1951.0 

2103.0 

2219.0 

1261.0 

SD 

46.1 

151.5 

27.8 

40.6 

105.2 

108.8 

89.8 

66.1 

CV(%) 

4.5 

4.3 

6.7 

5.8 

5.3 

5.2 

4.1 

5.3 

70%  PMT 

Mean 

13035.1 

44865.2 

5102.4 

8514.1 

25737.5 

26482.1 

29591.7 

15903.5 

Median 

12969.5 

45105.5 

5135.5 

8381.0 

28009.5 

28677.0 

29519.5 

15833,5 

SD 

537.4 

2384.0 

309.9 

671.6 

1550.0 

2254.0 

1514.2 

554.9 

CV(%) 

4.1 

5.4 

6.1 

7.9 

6.1 

8.5 

5.1 

3.5 

75%  PMT-Cy5,  70%  PMT-Cy3 
Mean  19371.4 

40990.5 

7277.1 

7487.8 

36491.5 

24824.8 

41453.4 

15200.2 

Median 

19359.5 

41011.5 

7366.0 

7491.0 

36285.5 

25081.0 

41727.0 

15295.0 

SD 

511.8 

2779.5 

485.2 

652.1 

1310.1 

1090.6 

2236.6 

954.6 

CV(%) 

2.6 

8.8 

8.7 

8.7 

3.6 

4.4 

5.4 

6.3 

Two  properties  regarding  the  RT-PCR  amplifica¬ 
tion  were  examined:  Do  the  amplified  probes  repre¬ 
sent  the  true  abundance  of  each  transcript  in  cells? 
Does  the  probe  saturate  the  arrays?  First,  we  mixed 
five  tk  genes  in  three  ratios  (1,  20,  and  400)  and  used 
them  as  PCR  templates.  Hybridizations  of  such  PCR 
product  to  a  small  five-tk  gene  array  showed  a  “com¬ 
pressed”  ratio  of  three-  to  fourfold  despite  a  400-fold 
input  template  ratio.  Nevertheless,  the  relative  order 
of  abundance  of  the  five  tk  templates  was  conserved. 
Saturation  has  likely  occurred  in  PCR,  which  may 
have  also  depleted  some  tk  primers  due  to  their  high 
degeneracy.  Second,  we  compared  single-probe  hy¬ 
bridizations  at  various  concentrations.  Four  to  six 
hundred  ng  of  PCR  amplified  products  typically  were 
used  for  random  prime  labeling.  The  labeled  probe 
was  purified  and  different  dilutions  were  applied  to 
the  arrays.  It  was  determined  that  one  tenth  of  the  la¬ 
beled  probe  per  channel  provided  optimal  signals  for 
analysis  without  saturating  the  arrays. 

Repeated  hybridizations  have  identified  distinctive 
patterns  among  thyroid,  breast,  and  prostate  cancer  cell 
lines  (Figure  1).  Although  the  panel  of  58  tk  genes  may 
not  represent  all  the  tk  expressed  in  any  single  cancer 


Figure  1  Multiple  hybridizations  confirmed  distinctive  patterns 
for  three  tumor  cell  lines.  (Array  1)  Prostate  cancer  cel!  lines 
Cy5:BPH-1-CAFTD-04/Cy3:BPH-1-CAFTD-02.  (Array  2)  Thyroid  cancer 
cell  lines  Cy5:FTC133/Cy3:FTC236.  (Array  3)  Breast  cancer  cel!  line 
Cy5:MCF7/Cy3:58  tk  mix. 


tissue,  we  have  reason  to  believe  they  represent  over 
50%  of  human  tk  genes  expressed  in  a  particular  tissue 
type  (Robinson  et  al.  1996).  These  distinctive  expres¬ 
sion  patterns,  as  unveiled  by  microarray,  become  molec¬ 
ular  fingerprints  of  tk  gene  expression  for  each  tumor. 

Acknowledgments 

Supported  by  a  grant  from  the  Director,  Office  of  Sci¬ 
ence,  Office  of  Biological  and  Environmental  Research  of 
the  US  Department  of  Energy  under  Contract  DE-AC03- 
76SF00098,  by  a  fellowship  from  the  Cancer  Research 
Foundation  of  America  (to  HBH),  and  by  grants  from  the 
Breast  Cancer  Research  Programs,  US  Army  Medical  Re¬ 
search  and  Material  Command,  US  Department  of  Defense 
(BC98-0937,  DC991395).  All  experiments  involving  human 
cells  or  cell  lines  were  approved  by  the  UC  Berkeley  and  the 
UCSF  Human  Subject  Use  Institutional  Review  Board. 


Literature  Cited 

Celis  JE,  Kruhoffer  M,  Gromova  I,  Frederiksen  C,  Ostergaard  M, 
Thykjaer  T,  Gromov  P,  Yu  J,  Palsdottir  H,  Magnusson  N,  Orn- 
toft  TF  (2000)  Gene  expression  profiling:  monitoring  transcrip¬ 
tion  and  translation  products  using  DNA  microarrays  and  pro- 
teomics.  FEBS  Lett  480:2-16 

Cheung  V,  Morley  M,  Aguilar  F,  Massami  A,  Kucherlapati  R, 
Childs  G  (1999)  Making  and  reading  microarrays.  Nature  Genet 
Suppl  21:15-19 

DeRisi  J,  Penland  L,  Brown  PO,  Bittner  ML,  Meltzer  PS,  Ray  M, 
Chen  Y,  Su  YA,  Trent  JM  (1996)  Use  of  a  cDNA  microarray  to 
analyze  gene  expression  patterns  in  human  cancer.  Nature  Genet 
14:457-460 

Iyer  VR,  Eisen  MB,  Ross  DT,  Schuler  G,  Moore  T,  Lee  JC,  Trent 
JM,  Staudt  LM,  Hudson  J,  Boguski  MS,  Lashkari  D,  Shalon  D, 
Botstein  D,  Brown  PO  (1999)  The  transcriptional  program  in  the 
response  of  human  fibroblasts  to  serum.  Science  283:83-86 

Robinson  D,  He  F,  Pretlow  T,  Kung  HJ  (1996)  A  tyrosine  kinase  pro¬ 
file  of  prostate  carcinoma.  Proc  Natl  Acad  Sci  USA  93:5958-5962 

Schena  M  (1996)  Genome  analysis  with  gene  expression  microar¬ 
rays.  BioEssays  18:427-431 


Appendix  2 


_ Volume  49(7):  925-926,  2001 

The  Journal  of  Histochemistry  &  Cytochemistry 

http://www.jhc.org 


BRIEF  REPORT 

Monitoring  Signal  Transduction  in  Cancer:  From  Chips  to  FISH 


Robert  A.  Lersch,  Jingly  Fung,  H.-Ben  Hsieh,  Jan  Smida,  and  Heinz-Ulrich  G.  Weier 

Life  Science  Division,  Lawrence  Berkeley  National  Laboratory,  Berkeley,  California  (RAL,H-BH,H-UGW);  School  of  Medicine, 
University  of  California,  San  Francisco,  California  (JF);  and  Arbeitsgruppe  Zytogenetik,  gsf-Forschungszentrum  fur  Umwelt 
und  Gesundheit,  Neuherberg,  Germany  (JS) 


SUMMARY  The  microarray  format  of  RNA  transcript  analysis  should  provide  new  clues  to 
carcinogenic  processes.  Because  of  the  complex  and  heterogeneous  nature  of  most  tumor 
samples,  histochemical  techniques,  particularly  RNA  fluorescent  in  situ  hybridization  (FISH), 
are  required  to  test  the  predictions  from  microarray  expression  experiments.  Here  we  de¬ 
scribe  our  approach  to  verify  new  microarray  data  by  examining  RNA  expression  levels  of  five 
to  seven  different  transcripts  in  a  very  few  cells  via  FISH.  (J  Histochem  Cytochem  49:925-926, 2001) 
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Our  laboratory  is  adapting  microarray  technology 
to  identify  new  genes  involved  in  the  formation  of  tu¬ 
mors  in  the  thyroid,  breast,  and  prostate.  Previously, 
the  transcription  levels  of  only  a  few  genes  could  be 
assayed  per  experiment.  Microarrays  circumvent  this 
limitation.  However,  a  second  problem  limits  progress 
in  cancer  research.  Most  tumors  are  a  mix  of  cell 
types,  due  in  part  to  the  normal  complexity  of  the  tis¬ 
sue  and  in  part  to  the  complexity  of  tumors  as  they 
evolve  from  benign  to  metastatic.  If  researchers  collect 
microarray  data  without  confronting  the  problem  of 
tumor  heterogeneity,  important  correlations  will  be 
missed. 

Preliminary  studies  performed  in  our  lab  and  else¬ 
where  indicated  that  solid  tumors  are  heterogeneous 
with  respect  to  oncogene  expression.  Figure  1  shows 
the  expression  of  brk  (O’Bryan  et  al.  1991)  in  a  child¬ 
hood  thyroid  cancer  that  arose  after  the  Chernobyl 
nuclear  accident  (Zitzelsberger  et  al.  1999).  In  this 
case,  two  cells  strongly  express  brk  (arrows  point  to 
their  nuclei)  and  a  nearby  cell  does  not  express  brk. 
These  studies  were  performed  with  touchprep  papil¬ 
lary  thyroid  tumor  specimens  using  filter-based  micro¬ 
scope  systems  and  only  two  different  cDNA  probes.  If 
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more  than  two  probes  are  used  per  experiment,  more 
slides,  time,  and  reagents  must  also  be  used.  All  exper¬ 
iments  involving  human  cells  or  cell  lines  were  ap¬ 
proved  by  the  U.C.  Berkeley  Human  Subject  IRB. 

To  address  this  problem,  we  are  developing  a  sys¬ 
tem  to  simultaneously  measure  cell-by-cell  RNA  levels 
of  several  different  markers.  The  proposed  scheme  to 
discriminate  benign  and  malignant  neoplasms  and  to 
identify  new  prognostic  markers  will  take  advantage 
of  new  methods  and  our  experience  with  oncogene  ac¬ 
tivation  in  thyroid  tumors.  Our  short-term  goal  is  to 
develop  a  system  which  determines  the  relative  level  of 
expression  of  five  tyrosine  kinase  genes  using  FISH- 
based  methods  and  Spectral  Imaging  (Sim).  Existing 
Sim  instrumentation  can  record  fluorescence  spectra 
from  400  nm  to  1100  nm  with  about  10-nm  resolu¬ 
tion,  whereas  the  resolution  of  a  light  microscope  is 
about  1  p.m. 

Sim  combines  the  techniques  of  fluorescence  mi¬ 
croscopy,  charge-coupled  device  (CCD)  camera,  and 
Fourier  spectroscopy.  The  light  emitted  from  each 
point  of  the  sample  is  collected  with  the  microscope 
objective  and  sent  to  a  collimating  lens.  This  light 
travels  through  an  optical  head  (interferometer),  is  fo¬ 
cused  on  a  CCD,  camera,  and  the  resulting  data  are 
processed  with  a  computer.  The  interferometer  divides 
each  incoming  beam  (light  from  the  microscope)  into 
two  coherent  beams  and  creates  a  variable  optical 
path  difference  (OPD)  between  them.  The  beams  are 
then  combined  to  interfere  with  each  other.  The  re¬ 
sulting  interference  intensity  is  measured  by  the  CCD 
detector  as  a  function  of  the  OPD.  The  intensity  vs 
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Figure  1  Frozen  tissue  touchprep  of  childhood  thyroid  tumor 
(S246)  was  hybridized  with  a  cDNA  probe  for  brk.  The  probe  was 
labeled  with  digoxigenin  and  detected  with  rhodamine-conju- 
gated  antibodies.  About  5%  of  the  cells  expressed  brk  at  high  lev¬ 
els,  the  other  cells  at  low  levels.  Arrows  indicate  two  high  brk- 
expressing  cells. 

OPD  is  an  “interferogram.”  The  spectrum,  i.e.,  inten¬ 
sity,  as  a  function  of  wavelength  can  be  recovered 
from  the  interferogram  using  Fourier  transforms.  The 
first  published  application  of  Sim  (spectral  karyo¬ 
typing  or  SKY)  (Schroeck  et  al.  1996)  was  developed 
to  screen  metaphase  spreads  for  translocations  with 
24  chromosome-specific  whole  chromosome  painting 
probes  labeled  with  Spectrum  Green,  Spectrum  Or¬ 
ange,  Texas  Red,  Cy5,  or  Cy5.5,  or  combinations 
thereof.  The  multiple  bandpass  filter  set  (Chroma 
Technology;  Brattleboro,  VT)  used  for  fluorochrome 
excitation  was  designed  to  provide  broad  emission 
bands  (giving  a  fractional  spectral  reading  from  ~450 
nm  to  ~850  nm). 

Although  we  can  use  much  of  the  existing  technol¬ 
ogy,  the  ratio-labeling  color  scheme  for  SKY  will  not 
work  for  RNA  detection  because  different  labeled 
cDNA  probes  might  co-localize.  Therefore,  we  will  la¬ 
bel  each  cDNA  probe  with  a  unique  reporter.  The  flu¬ 
orescence  spectra  of  the  different  reporter  molecules 
can  partially  overlap,  because  the  signals  can  be  re¬ 
solved  by  a  computer  algorithm  termed  “Spectral  Un- 
Mixing  (SUN)”  (Applied  Spectral  Imaging;  Carlsbad, 
CA).  SUN  enables  us  to  deconvolute  overlapping  spec¬ 
tra  and  recover  single  component  images  from  the 
spectral  image. 

The  existing  commercial  SKY  metaphase  chromo¬ 
some  analysis  software  will  be  modified  to  increase  its 
automated  signal  processing,  RNA  identification  in  in¬ 
terphase  cells,  integration  of  cDNA  probe  signals,  and 
databasing  of  results.  We  will  only  refine  the  analysis 
software;  the  acquisition  methods  of  capturing  a  spec¬ 
tral  image  and  a  separate  high-contrast  monochrome 
DAPI  image  will  remain  unchanged.  A  prototype  of 


this  system  for  analysis  of  protein  markers  has  already 
been  developed  (Tsurui  et  al.  2000)  and  will  aid  our 
efforts  to  analyze  RNA  transcript  levels. 

Initially  we  will  hybridize  cDNA  probes  to  brk,  ret, 
c-met,  trk,  axl/ufo  (obtained  from  Research  Genetics; 
Huntsville,  AL)  using  five  cyanine  dyes  (Amersham; 
Arlington  Heights,  IL).  On  the  basis  of  previous  mi¬ 
croarray  results,  we  will  also  add  members  of  the  Eph 
tyrosine  kinase  family.  cDNA  probes  will  be  prepared 
by  incorporating  fluorochrome-labeled  deoxynucleo- 
side  triphosphates  by  random  priming  or  PCR  ampli¬ 
fication  (Weier  et  al.  1994,19 95).  We  will  also  build  a 
reference  spectra  database  for  these  probes.  Standard 
fixation  and  hybridization  protocols  should  detect  the 
target  RNAs.  It  is  unclear  if  we  will  need  to  prepare 
custom  blocking  agents  to  DNA  repeats  in  the  3'- 
untranslated  region  of  some  of  our  targets.  Using  arti¬ 
ficial  mixtures  of  existing  cell  lines,  we  will  develop 
software  modules  needed  to  measure  intracellular  lev¬ 
els  of  five  RNA  species  and  to  determine  important 
process  parameters:  sensitivity,  accuracy,  and  repro¬ 
ducibility. 
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SUMMARY  Abnormal  expression  of  tyrosine  kinase  (TK)  genes  is  common  in  tumors,  in 
which  it  is  believed  to  alter  cell  growth  and  response  to  external  stimuli  such  as  growth  fac¬ 
tors  and  hormones.  Although  the  etiology  and  pathogenesis  of  carcinomas  of  the  thyroid 
or  breast  remain  unclear,  there  is  evidence  that  the  expression  of  TK  genes,  such  as  recep¬ 
tor  tyrosine  kinases,  or  mitogen-activated  protein  kinases,  is  dysregulated  in  these  tumors, 
and  that  overexpression  of  particular  TK  genes  due  to  gene  amplification,  changes  in  gene 
regulation,  or  structural  alterations  leads  to  oncogenic  transformation  of  epithelial  cells. 
We  developed  a  rapid  scheme  to  measure  semiquantitatively  the  expression  levels  of  50- 
100  TK  genes.  Our  assay  is  based  on  RT-PCR  with  mixed  based  primers  that  anneal  to  con¬ 
served  regions  in  the  catalytic  domain  of  TK  genes  to  generate  gene-specific  fragments. 
PCR  products  are  then  labeled  by  random  priming  and  hybridized  to  DNA  microarrays  car¬ 
rying  known  TK  gene  targets.  Inclusion  of  differently  labeled  fragments  from  reference  or 
normal  cells  allows  identification  of  TK  genes  that  show  altered  expression  levels  during 
malignant  transformation  or  tumor  progression.  Examples  demonstrate  how  this  innova¬ 
tive  assay  might  help  to  define  new  markers  for  tumor  progression  and  potential  targets 
for  disease  intervention,  (j  Histochem  Cytochem  49:673-674, 2001) 
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Because  the  malignant  transformation  of  epithelial 
cells  and  progression  of  carcinomas  are  accompanied 
by  changes  in  the  expression  of  receptor  and  cytosolic 
tyrosine  kinase  (TK)  genes,  we  set  out  to  develop  an 
innovative  DNA  microarray-based  assay  to  simulta¬ 
neously  determine  the  relative  expression  level  of  50- 
100  TK  genes  using  just  a  few  cells. 

A  number  of  studies  have  shown  that  tumor  develop¬ 
ment  is  accompanied  by  at  least  two  changes:  (a)  a 
change  in  the  way  cells  interact  with  their  environment 
via  membrane-bound  receptors,  and  (b)  a  change  in  how 
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signals  originating  from  these  receptors  are  transduced 
from  the  cell  membrane  to  the  cytoplasm  and  the  nu¬ 
cleus.  Among  the  hundreds  of  genes  involved  in  recep¬ 
tor-mediated  signal  transduction,  only  a  few  are  aber¬ 
rantly  expressed  in  tumors.  A  major  motif  in  signal 
transduction  is  the  selective  phosphorylation  and  de¬ 
phosphorylation  of  tyrosine  residues  in  protein  factors 
involved  in  signal  processing.  Proteins  that  phosphory- 
late  tyrosine  residues  are  products  of  genes  belonging  to 
the  family  of  TK  genes.  The  number  of  known  TK  genes 
has  grown  steadily  in  recent  years,  and  the  temporal  and 
tissue-specific  expression  of  ~  100  different  TK  genes  in 
normal  cells  is  carefully  orchestrated.  The  level  of  ex¬ 
pression  of  certain  TK  genes  is  increased  in  many  human 
tumors  (Luttrell  et  al.  1994;  Liu  et  al.  1995).  Simple 
overexpression  of  TK  genes  due  to  gene  amplification  or 
to  changes  in  the  regulation  of  gene  expression  may  lead 
to  oncogenic  transformation.  This  has  been  clearly  doc¬ 
umented  for  the  erbB2  protein,  the  product  of  the  Her- 
2/neu  proto-oncogene  and  other  members  of  the  erb  B 
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family  (Tse  et  al.  1997).  In  addition,  many  tumors  have 
acquired  structurally  altered  TK  proteins  or  abnormal 
expression  patterns  through  de  novo  mutational  events. 
When  chromosomes  became  rearranged,  the  catalytic 
domain  of  a  TK  gene  was  found  fused  to  the  amino 
terminal  of  another  protein,  thus  creating  a  new  trans¬ 
forming  activity  as  well  as  a  new  expression  pattern. 
Well-known  examples  of  this  mechanism  of  oncogene 
activation  are  the  bcr/abi-fusion  protein  in  chronic  mye¬ 
loid  leukemia  with  a  translocation,  t(9;22),  and  the  acti¬ 
vation  of  the  receptor  TKs  ret  and  trk  in  papillary  thyroid 
cancer  (Sozzi  et  al.  1992;  Jhiang  and  Mazzaferri  1994). 

Various  protein  factors  can  be  mis-expressed  and,  in 
combination  with  other  events,  might  constitute  one  of 
several  factors  leading  to  the  onset  and/or  progression 
of  cancer.  Factors  including  cell  cycle-specific  enzymes, 
hormone  receptors,  and  peptide  growth  factors  have 
been  reported  as  having  prognostic  significance  in 
some  cases  of  prostate  cancer.  Overexpression  of  par¬ 
ticular  receptor  TK  genes  such  as  the  insulin-like 
growth  factor  receptors  (IGF-IRs),  the  epidermal 
growth  factor  receptor  (EGFR  or  erb  B)  family  of  re¬ 
ceptors,  focal  adhesion  kinase  (FAK),  or  the  proto¬ 
oncogenes  ret  and  Nyk/mer  have  been  shown  to  corre¬ 
late  with  progression  to  a  more  malignant  phenotype  in 
a  variety  of  tumors  (Resnik  et  al.  1998),  among  them 
carcinoma  of  the  prostate  (Ling  and  Kung  1995;  Daw¬ 
son  et  al.  1998).  Detailed  knowledge  about  TK  gene  ex¬ 
pression  and  its  relation  to  tumor  progression  might  in¬ 
crease  our  understanding  of  how  tumors  grow  and  help 
us  design  assays  to  more  accurately  stage  tumors. 

Our  project  targets  the  development  of  a  novel  as¬ 
say  format  that  enables  us  to  determine  the  level  of  ex¬ 
pression  of  many  different  genes.  A  rapid  assay  uses 
DNA  microarrays  carrying  small  amounts  of  individ¬ 
ual  TK  gene-specific  targets  to  simultaneously  deter¬ 
mine  the  expression  level  of  up  to  100  TK  genes  using 
a  small  number  of  cells.  We  cloned  and  characterized 
TK  genes  expressed  in  thyroid  cancers  and  in  seven  dif¬ 
ferent  breast  cancer  cell  lines.  Using  mixed-base  oligo¬ 
nucleotides  specific  for  conserved  domains  in  the  cata¬ 
lytic  domains  of  TK  genes,  our  PCR  assays  amplified 
^15 9—1 71-bp  fragments  of  expressed  TK  genes.  We 
size-selected  and  cloned  the  PCR  products  into  plas¬ 
mids.  As  of  August  2000,  we  had  identified  more  than 
50  TK  genes  expressed  in  thyroid  cancers  plus  an  addi¬ 
tional  172  TK  fragment  clones  derived  from  breast 
cancer  cell  lines  and  ~250  TK  fragment  clones  from 
radiation-induced  thyroid  cancers  in  the  prescreening 
and  sequencing  steps.  Previously,  we  finished  the  con¬ 
struction  of  a  robotic  system  to  print  DNA  microar¬ 
rays  with  about  100  sequences  on  glass  slides  (unpub¬ 
lished  data).  The  performance  of  the  system  was  tested 
by  hybridization  of  fluorochrome-labeled  TK  gene-spe¬ 
cific  PCR  fragments  onto  the  TK-specific  DNA  mi¬ 
croarrays.  All  experiments  involving  human  cells  or 


cell  lines  were  approved  by  the  U.C.  Berkeley  Human 
Subject  Use  Institutional  Review  Board.  These  experi¬ 
ments  enabled  us  to  optimize  hybridization  and  wash¬ 
ing  conditipns  and  to  generate  data  regarding  the  rela¬ 
tive  level  of  gene  expression.  For  example,  the  184A1 
and  184A1TH  cell  lines  are  closely  related  non-tumor- 
igenic  and  tumorigenic  human  mammary  epithelial  cell 
lines  derived  from  the  same  normal  breast  tissue  speci¬ 
men  and  transformed  in  vitro.  After  hybridizing  Cy5- 
labeled  TK  fragments  prepared  from  the  cell  line 
184A1  and  Cy3-labeled  TK  fragments  prepared  from 
cell  line  184A1TH  to  a  DNA  microarray  carrying  58 
TK  gene  fragments  and  displaying  Cy3  fluorescence 
signals  in  green  and  Cy5  signals  in  red,  differences  be¬ 
tween  the  cell  lines  became  readily  apparent.  Genes  ex¬ 
pressed  at  a  higher  level  in  184A1TH  cells  led  to  in¬ 
creased  green  signals  on  the  array,  while  those  genes 
whose  expression  level  is  lower  in  184A1TH  cells 
compared  to  184A1  cells  generated  spots  that  exhib¬ 
ited  stronger  red  fluorescence.  Similarly,  differences  be¬ 
tween  thyroid  tumor  cell  lines  could  be  demonstrated 
in  simple  dual-color  hybridization  experiments. 
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Abstract 

Numerical  chromosome  aberrations  are  incompatible  with  normal  human  development.  Our  laboratories  develop  hybridization- 
based  screening  tools  that  generate  a  maximum  of  cytogenetic  information  for  each  polar  body  or  blastomere  analyzed.  The 
methods  are  developed  considering  that  the  abnormality  might  require  preparation  of  case-specific  probes  and  that  only  one  or 
two  cells  will  be  available  for  diagnosis,  most  of  which  might  be  in  the  interphase  stage.  Furthermore,  assay  efficiencies  have  to 
be  high,  since  there  is  typically  not  enough  time  to  repeat  an  experiment  or  reconfirm  a  result  prior  to  fertilization  or  embryo 
transfer.  Structural  alterations  are  delineated  with  breakpoint-spanning  probes.  When  screening  for  numerical  abnormalities,  we 
apply  a  Spectral  Imaging-based  approach  to  simultaneously  score  as  many  as  ten  different  chromosome  types  in  individual 
interphase  cells.  Finally,  DNA  micro-arrays  are  under  development  to  score  all  of  the  human  chromosomes  in  a  single  experiment 
and  to  increase  the  resolution  with  which  micro-deletions  can  be  delineated.  ©  2001  Elsevier  Science  Ireland  Ltd.  All  rights 
reserved. 

Keywords:  Cytogenetics;  Chromosomes;  Hybridization;  FISH;  DNA  micro-arrays;  Digital  image  analysis  (blastomeres) 


t.  Introduction 

Carriers  of  balanced  translocations  have  an  elevated 
risk  of  producing  aneuploid  germ  cells  due  to  disturbed 
homologue  pairing.  The  resulting  partial  or  total  aneu- 
somies  lead  to  spontaneous  abortions,  stillbirth  or 
severe  deficiencies  and  disease.  Assisted  reproductive 
technology  now  offers  couples  at  risk  several  diagnostic 
approaches  to  reduce  the  risk  of  carrying  an  affected 
fetus.  If  the  woman  carries  the  abnormality,  first  polar 
bodies  can  be  analyzed  immediately  after  oocyte  har¬ 
vest.  Following  in  vitro  fertilization,  pre-implantation 
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genetic  analysis  (PGD)  can  be  performed  on  individual 
blastomeres  biopsied  from  3-day-old  embryos.  Since 
most  of  the  embryonic  cells  will  be  found  in  interphase 
stage,  the  diagnostic  approach  will  have  to  work  reli¬ 
ably  with  either  the  less  condensed  chromatin  in  inter¬ 
phase  cell  nuclei  or  the  highly  condensed  DNA  in  polar 
bodies  (PB’s). 

Our  collaborating  laboratories  have  long  been  in¬ 
volved  in  the  development  of  nucleic  acid  hybridiza¬ 
tion-based  procedures  for  the  rapid  detection  of 
structural  and  numerical  chromosome  abnormalities. 
Here,  we  report  the  present  state  of  hybridization-based 
technologies  for  interphase  cell  analysis  in  PGD.  Since 
only  one  or  two  ceils  are  available  for  analysis,  our 
approaches  are  geared  towards  obtaining  a  maximum 
of  cytogenetic  information  per  experiment. 
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2.  Material  and  methods 

„  For  more  than  a  decade,  our  laboratories  have  been 
involved  in  the  development  of  technologies  for  analy¬ 
sis  of  interphase  and  metaphase  cells.  In  collaboration 
with  scientists  at  the  University  of  California,  San 
Francisco,  and  the  St  Barnabas  Medical  Center,  Liv¬ 
ingston,  researchers  at  the  E.O.  Lawrence  Berkeley 
National  Laboratory  study  the  chromosomal  composi¬ 
tion  of  blastomeres  with  regard  to  numerical  as  well  as 
structural  aberrations.  The  technical  aspects  of  our 
probe  preparation  and  multicolor  detection  protocols 
have  been  published  previously  (Weier  et  al.,  1994; 
Jossart  et  al.,  1996;  Cassel  et  al.,  1997). 

A  major  goal  of  our  technical  developments  is  to 
maximize  the  number  of  chromosomal  targets  that  can 
be  scored  simultaneously.  Briefly,  probes  specific  for 
repeated  DNA  on  chromosomes  15,  X,  and  Y  pur¬ 
chased  from  Vysis  (Downers  Grove,  IL)  were  labeled 
with  either  a  green  or  red  fluorochrome  (Spectrum 
Green  or  Spectrum  Orange,  respectively).  The  probes 
specific  for  chromosome  9,13,14,16,18,21,  and  22 
were  prepared  in  house  and  labeled  by  random  priming 
(BioPrime  Kit,  GIBCO/LTI,  Gaithersburg,  MD)  incor¬ 
porating  biotin- 14-dCTP  (part  of  the  BioPrime  Kit), 
digoxigcnin-1 1-dUTP  (Roche  Molecular  Biochemicals, 
Indianapolis,  IN),  fluorescein- 12-dUTP  (Roche  Molec¬ 
ular  Biochcmicals)  (Weier  et  al.,  1995),  or  Cy3-dUTP 
(Amcrsham,  Arlington  Heights,  IN).  Bound  biotiny¬ 
lated  probes  were  detected  with  avidin-Cy5,  and  bound 
digoxigcnin-labcled  probes  were  detected  with  Cy5.5- 
conjugated  antibodies  against  digoxin  (Sigma,  St. 
Louis,  MO).  Between  0.5  and  3  pi  of  each  probe  along 
with  1  pi  human  COT1™  DNA  (1  mg/ml,  GIBCO/ 
LTI)  and  1  pi  salmon  sperm  DNA  (20  mg/ml,  3'-5\ 
Boulder,  CO)  were  precipitated  with  1  pi  glycogen 
(Roche  Molecular  Biochemicals,  1  mg/ml)  and  1/10 
volume  of  3  M  sodium  acetate  in  2  volumes  of  2- 
propanol,  air  dried  and  resuspended  in  3  pi  water, 
before  7  pi  of  hybridization  master  mix  [78.6%  for- 
mamidc  (FA,  GIBCO/LTI),  14.3%  dextran  sulfate  in 
2.9  x  SSC,  pH  7.0  (1  x  SSC  is  150  mM  NaCl,  15  mM 
Na  citrate)]  were  added.  This  gave  a  total  hybridization 
mixture  of  10  pi. 

All  blastomeres  used  in  the  probe  developments  were 
obtained  from  embryos  donated  by  patients  enrolled  in 
the  IVF  Programs  of  The  University  of  California,  San 
Francisco,  or  The  Institute  for  Reproductive  Medicine 
and  Science  of  Saint  Barnabas  Medical  Center.  In 
accordance  with  guidelines  set  by  the  internal  review 
boards  of  these  Medical  Centers,  written  consent  was 
obtained  from  the  patients  in  each  case.  Embryo  biop¬ 
sies  and  blastomere  fixations  were  carried  out  as  de¬ 
scribed  (Munne  et  al.,  1994,  1996).  As  indicated  below, 
embryos  used  for  some  studies  had  arrested  develop¬ 
ment  or  were  morphologically  abnormal. 


3.  Results  and  discussion 

Traditional  filter  based  microscope  systems  limited 
fluorescence  in  situ  hybridization  (FISH)  experiments 
to  the  simultaneous  use  of  typically  no  more  than  three 
to  five  differently  labeled  probes  for  interphase  analysis 
(Munne  and  Weier,  1996).  This  is  sufficient  to  detect 
structural  alterations  in  interphase  and  metaphase  cells 
or  score  a  few  chromosomes  in  interphase  cells  (Munne 
et  al.,  1994;  Munne  and  Weier,  1996;  Munne  et  al., 

1996) .  We  prepared  case-specific  breakpoint-spanning 
probe  contigs  to  identify  intrachromosomal  rearrange¬ 
ments  such  as  inversions  or  deletions  (Cassel  et  al., 

1997) .  The  same  approach  can  be  used  to  demonstrate 
interchromosomal  rearrangements  such  as  reciprocal 
translocations  (Munne  et  al.,  1998a;  Fung  et  al.,  1999; 
Weier  et  al.,  1999).  The  case-specific  probes  allow  one 
to  discriminate  between  a  normal  karyotype,  aneuploid 
cells,  and  a  balanced  karyotype  carrying  the  derivative 
chromosomes.  A  less  time  consuming  and,  thus,  less 
expensive  approach  using  DNA  probes  that  bind  distal 
to  the  respective  breakpoints  can  only  be  used  to  count 
the  number  of  chromosome  copies  and  thus  cannot 
discriminate  between  the  normal  and  the  balanced 
karyotypes  (Munne  et  al.,  1998b). 

3. 1 .  Detection  of  structural  chromosome  aberrations 

Our  scheme  for  the  detection  of  structural  alteration 
is  based  on  the  preparation  and  hybridization  of  two 
differently  labeled  DNA  probes  that  bind  on  both  sides 
of  the  respective  chromosome  breakpoints  (Cassel  et 
al.,  1997).  One  probe  will  be  detected  in  the  green 
fluorescence  wavelength  interval,  while  the  second 
probe  is  made  such  that  it  fluoresces  red.  Normal 
homologues  lack  the  rearrangement  and  the  probes 
produce  large  hybridization  domains  that  appear  either 
red  or  green  in  the  fluorescence  microscope.  At  the 
same  time,  we  counterstain  the  DNA  with  4,6-diamino- 
2-phenylindole  (DAPI)  which  fluoresces  blue  under  ul¬ 
traviolet  light  excitation.  If  the  cell  contains  a  derivative 
chromosome,  the  hybridization  result  will  show  a  red/ 
green  associated  or  partially  overlapping  signal  indica¬ 
tive  of  the  translocation  event.  Thus,  we  detect 
structural  alteration  and  score  homologues  at  the  same 
time. 

Our  work  is  greatly  facilitated  by  access  to  resources 
created  in  the  course  of  the  International  Human 
Genome  Sequencing  Project,  such  as  large  insert  ge¬ 
nomic  DNA  libraries  (bacterial  or  yeast  artificial  chro¬ 
mosomes  (BACs  or  YACs,  respectively),  high 
resolution  physical  maps  or  collections  of  cytogeneti¬ 
cally  mapped  DNA  probes  (Chen  et  al.,  1996;  Kim  et 
al.,  1996;  Korenberg  et  al.,  1999).  This  enables  us  to 
prepare  case-specific  probe  sets  suitable  for  interphase 
cell  analyses  of  most  patient  cells  within  a  few  weeks. 
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Once  optimized,  these  probe  sets  allow  to  rapidly  deter¬ 
mine  the  exact  number  of  normal  chromosomes  and 
^derivative  chromosomes  in  somatic  cells  from  translo¬ 
cation  carriers  as  well  as  their  germ  cells  or  offspring. 
So  far,  however,  these  procedures  failed  to  produce  the 
desired  increase  in  pregnancy  rates  in  cases  where  one 
spouse  carried  a  balanced  reciprocal  translocation.  Our 
concern  is  that  the  impaired  homologue  pairing  in  the 
carriers  leads  to  gain  or  loss  of  other  chromosomes 
which  remains  undetected  in  assays  scoring  only  the 
translocation  chromosomes. 

The  recent  introduction  of  Spectral  Imaging  (Sim) 
now  allows  one  to  interrogate  a  much  larger  number  of 
targets,  thus  producing  a  more  comprehensive  picture 
of  the  chromosomal  composition  of  the  cells.  Sim 
allows  an  investigator  to  discriminate  an  increased 
number  of  fluorescent  probes  by  exciting  fluorescent 
molecules  over  a  broad  spectral  range  and  by  recording 
the  fluorescence  emission  spectral  using  an 
interferometer. 

3.2.  Detection  of  numerical  chromosome  abnormalities 
using  Spectral  Imaging  (Sim) 

Chromosome  abnormalities  occur  with  astonishing 
frequency  in  humans,  being  present  in  an  estimated 
10-30%  of  all  fertilized  eggs.  Over  25%  of  all  the 
miscarriages  are  monosomic  or  trisomics,  making  aneu- 
ploidy  the  leading  known  cause  of  pregnancy  loss. 
Ideally,  one  likes  to  detect  ancuploidy  involving  any  of 
the  24  human  chromosomes  for  preimplanation  genetic 
and  prenatal  diagnosis.  Thus,  an  analytical  method  to 
enumerate  as  many  chromosomes  as  possible  in  few 
interphase  cells  is  highly  desirable.  Using  a  set  com¬ 
prized  of  seven  chromosome-specific  probes  (chromo¬ 
some  10,  14,  16,  18,  22,  X  and  Y)  hybridized  to 
lymphocyte  interphase  nuclei,  we  demonstrated  that 
Spectral  Imaging  system  provides  a  significant  improve¬ 
ment  over  conventional  filter-base  microscope  systems 
for  enumeration  of  multiple  chromosomes  in  interphase 
nuclei  (Fung  et  al.,  1998b). 

Using  mostly  yeast  or  bacterial  artificial  chromosome 
probes  for  cytogenetic  analyses  of  blastomeres  and 
detection  of  structural  alterations,  we  are  building  pan¬ 
els  of  probes  to  simultaneously  score  10  or  more  differ¬ 
ent  chromosomes.  Further  increases  in  the  number  of 
probes  is  complicated  due  to  occasional  overlap  of 
chromosome  domains  or  local  variation  in  hybridiza¬ 
tion  efficiency.  We  developed  a  10-chromosomes  probe 
set  (chromosomes  9,  13,  14,  15,  16,  18,  21,  22,  X  and  Y) 
for  the  purpose  of  labeling  DNA  targets  most  fre¬ 
quently  associated  with  aneuploidy  and  spontaneous 
abortions  and  tested  its  application  in  PGD  (blas¬ 
tomeres  from  abnormal  human  preimplantation  em¬ 
bryos)  and  prenatal  diagnosis  (uncultured  amniocytes 
obtained  by  amniocentesis).  Results  demonstrated  in¬ 


creasing  levels  of  background  fluorescence  on  different 
cells  after  hybridization  in  the  order:  (uncultured  am¬ 
niocytes)  » (blastomeres)  >  (interphase  cells  from 
lymphocytes).  All  blastomeres  fixed  for  this  study  {N  — 
25)  spread  very  well.  Fourteen  nuclei  (56%)  showed 
interpretable  hybridization  results,  and  most  of  them 
were  karyotyped  as  abnormal,  since  all  those  cells  were 
from  1  and  3  PN  human  embryos,  and  had  arrested 
development  or  were  morphologically  abnormal.  The 
signals  from  11  nuclei  (44%)  were  faint.  This  may  be 
related  to  the  quality  of  the  embryos,  since  all  of  them 
were  developing  abnormally.  The  fixation  of  uncultured 
amniocytes  on  slides  for  Spectral  Imaging  analysis 
turned  out  to  be  somewhat  difficult.  Most  nuclei  were 
not  flattened  out,  presenting  a  problem  due  to  the 
limited  focal  depth.  Overlapping  signal  domains  were  a 
problem  in  uncultured  amniocytes,  where  only  about 
20%  of  all  cells  showed  interpretable  spreads.  In  sum¬ 
mary,  Spectral  Imaging  has  demonstrated  advantages 
for  evaluating  numerical  chromosomal  abnormalies  in 
single  interphase  cells.  Its  utility  for  chromosome  scor¬ 
ing,  however,  remains  limited  due  to  chromosome  do¬ 
main  overlap. 

3.3.  DNA  micro -arrays  (chips) 

The  DNA  micro-arrays  represent  an  exciting  new 
technology  with  applications  ranging  from  gene  expres¬ 
sion  profiling  to  determination  of  gene  copy  number 
changes  in  tumors.  We  are  presently  investigating  this 
approach  in  which  the  DNA  probes  are  immobilized  on 
glass  slides  as  a  strategy  complementing  FISH  studies. 
The  DNA  from  the  embryonic  cells  is  labeled  in  one 
color  (e.g.  red),  while  an  equal  amount  of  a  reference 
DNA  probe  is  labeled  in  a  different  color  (e.g.  green). 
These  differentially  labeled  DNA  are  combined,  dena¬ 
tured  and  hybridized  to  a  DNA  micro-array  or  ‘chip’  in 
a  quantitative  manner.  Results  are  obtained  after  read¬ 
ing  the  micro-arrays  with  specially  designed  fluores¬ 
cence  scanners  as  red/green  or  red/infrared  fluorescence 
ratios.  After  normalization,  every  increase  or  decrease 
from  the  average  ratio  indicates  an  abnormal  number 
of  copies  of  the  hybridization  target. 

In  practice,  the  DNA  contained  in  a  single  cell  is  not 
sufficient  to  generate  measurable  signals.  The  com¬ 
monly  used  protocols  therefore  include  a  DNA  in  vitro 
amplification  step  using  a  random  primer  or  oligonucle¬ 
otide  with  arbitrary  sequence  prior  to  labelling.  The 
DNA  to  be  immobilized  can  be  obtained  by  standard 
isolation  protocols  or  by  in  vitro  DNA  amplification. 
We  use  a  DNA  spotter  based  on  the  design  published 
by  Brown’s  group  at  Stanford  University  (Schena  et  al., 
1996).  This  allows  us  to  spot  small  amounts  of  DNA  on 
poly-L-lysine  coated  glass  slides  with  a  100-200  micron 
pitch.  A  288-spot  DNA  micro-array  like  the  test  array 
depicted  schematically  in  Fig.  1  then  measures  no  more 
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Fig.  I.  Prototype  DNA  micro-array  for  technology  development.  The 
chip  contains  from  left  to  right:  control  DNA  spots  (six  samples), 
PCR  amplified  YAC  DNA  (12  clones),  BAC  array  1  (46  clones  for 
copy  number  determination  in  breast  cancer  research),  cDNA  (62 
tyrosine  kinase  cDNA  for  expression  profiling),  and  BAC  array  2  (24 
clones  for  chromosome  enumeration  in  PGD).  All  samples  are  spot¬ 
ted  in  duplicate.  The  total  number  of  DNA  spots  of  300  arrayed  with 
a  180  pm  pitch  leads  to  array  dimensions  of  2.08  x  4.42  mm2.  Each 
slide  contains  two  identical  micro-arrays. 

than  a  few  square  millimeters.  This  small  size  of  micro¬ 
arrays  is  an  advantage  over  larger  arrays  prepared  on 
nylon  filters,  because  it  requires  less  amount  of  labeled 
DNA  sample  in  the  hybridization  reaction. 

We  are  presently  preparing  DNA  test  chips  to  de¬ 
velop  the  methods  and  study  parameters  such  as  probe 
preparation,  hybridization  conditions  and  chip  reader 
performance.  Our  long  term  objective  is  the  develop¬ 
ment  of  reliable  procedures  to  detect  structural  as  well 
as  numerical  abnormalities  using  a  combination  of 
‘FISH  and  chip’  technology.  Chips  to  be  used  in  those 
studies  will  carry  several  hybridization  targets  per  chro¬ 
mosome  arm,  thus  allowing  a  more  detailed  gene 
dosage  determination  or  delineation  of  full  or  partial 
aneusomies  than  the  FISH-based  assays. 
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SUMMARY  In  vitro  model  cell  systems  are  important  tools  for  studying  mechanisms  of 
radiation-induced  neoplastic  transformation  of  human  epithelial  cells.  In  our  study,  the 
human  thyroid  epithelial  cell  line  HTori-3  was  analyzed  cytogenetically  following  expo¬ 
sure  to  different  doses  of  a-  and  y-irradiation  and  subsequent  tumor  formation  in  a  thymic 
nude  mice.  Combining  results  from  G-banding,  comparative  genomic  hybridization,  and 
spectral  karyotyping,  chromosome  abnormalities  could  be  depicted  in  the  parental  line 
HTori-3  and  in  nine  different  HTori  lines  established  from  the  developed  tumors.  A 
number  of  chromosomal  aberrations  w  ere  found  to  be  characteristic  for  simian  virus  40 
immortalization  and/or  radiation-induced  transformation  of  human  thyroid  epithelial 
cells.  Common  chromosomal  changes  in  cell  lines  originating  from  different  irradiation 
experiments  were  loss  of  8q23  and  13cen-q21  as  well  as  gain  of  lq32-qter  and  2qll.2- 
q  14.1.  By  comparison  of  chromosomal  aberrations  in  cell  lines  exhibiting  a  different  tu- 
morigenic  behavior,  cytogenetic  markers  important  for  the  tumorigenic  process  were 
studied.  It  appeared  that  deletions  on  chromosomes  9q32-q34  and  7q21-q31  as  well  as  an 
increased  copy  number  of  chromosome  20  were  important  for  the  tumorigenic  phenotype. 
A  comparative  breakpoint  analysis  of  the  marker  chromosomes  found  and  those  observed 
in  radiation-induced  childhood  thyroid  tumors  from  Belarus  revealed  a  coincidence  for  a 
number  of  chromosome  bands.  Thus,  the  data  support  the  usefulness  of  the  established 
cell  system  as  an  in  vitro  model  to  study  important  steps  during  radiation-induced  malig¬ 
nant  transformation  in  human  thyroid  cells.  ©  2001  Wiley-Liss,  Inc. 
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INTRODUCTION 

The  neoplastic  transformation  of  normal  human  ep¬ 
ithelial  cells  to  cancer  cells  is  a  multistep  process 
associated  with  immortalization,  aberrant  growth 
control,  and  malignancy  [1],  For  a  detailed  under¬ 
standing  of  this  multistage  nature  of  human  carci¬ 
nogenesis.  it  is  essential  to  develop  in  vitro  models 
of  human  epithelial  cell  transformation.  So  far.  a 
number  of  human  epithelial  cell  lines  from  skin, 
bronchus,  thyroid  gland,  prostate,  and  ureteral  tis¬ 
sue  have  been  immortalized  using  viral  constructs. 
In  four  of  these  cel)  lines  1 1—7].  tumor  formation  in 
athymic  nude  mice  was  observed  following  irradia¬ 
tion  in  vitro.  For  the  immortalization  process,  viral 
vectors  (simian  virus  40  [SV40].  Epstein  Barr,  ad¬ 
enovirus,  and  human  papillomavirus)  have  been 
used  that  function  through  the  dysregulation  of  the 
p53  gene  and  DNA  damage-induced  checkpoints. 
Cell  lines  immortalized  in  this  way  are  often  ge¬ 
netically  unstable  and  exhibit  complex  rearranged 
karyotypes  complicating  the  identification  of  cyto¬ 
genetic  markers  related  to  radiocarcinogenesis. 
G-banding  analysis  is  not  sufficient  to  resolve  com¬ 
plex  changes  and  also  to  detect  hidden  abnormali¬ 
ties  so,  therefore,  molecular  cytogenetic  techniques 
such  as  spectral  karyotyping  (SKY )  and  compara¬ 
tive  genomic  hybridization  (CGH)  must  be  addi¬ 
tionally  applied.  SKY  has  been  developed  to  char¬ 
acterize  human  metaphase  spreads  after 
simultaneously  hybridizing  the  chromosomes  with 
24  chromosome-specific  whole  chromosome  paint¬ 
ing  probes  [8].  This  combination  of  spectroscopic 
detection  methods  and  in  situ  hybridization  assays 
using  fluorescently  tagged  probes  allows  a  simul¬ 
taneous  depiction  and  analysis  of  the  whole  human 
chromosome  complement.  Interchromosomal  rear¬ 
rangements  can  be  recognized  as  differentially 
stained  chromosome  segments,  thus  greatly  facili¬ 
tating  the  characterization  of  structural  aberrations, 
especially  those  with  complex  patterns  (8—10]  and 
of  hidden  chromosome  abnormalities  [11,12],  If 
this  technique  is  used  in  addition  to  comparative 
genomic  hybridization  [13],  which  depicts  copy 
number  changes  with  high  accuracy,  a  comprehen¬ 
sive  picture  of  chromosomal  abnormalities  devel¬ 
oped  during  the  neoplastic  cell  transformation  can 
be  obtained. 

In  our  study,  SV40-immortaIized  and  radia¬ 
tion-transformed  human  thyroid  cell  lines  (HTori) 
were  analyzed  cytogenetically  by  G-banding.  SKY 


and  CGH.  We  investigated  the  immortalized  paren¬ 
tal  cell  line  as  well  as  nine  radiation-transformed 
(a-  and  y-irradiated  with  different  doses)  cell  lines 
and  identified  specific  chromosomal  changes  at¬ 
tributable  to  the  immortalization  process,  the  tu- 
morigenic  behavior,  and/or  the  different  radiation 
qualities. 

MATERIALS  AND  METHODS 
Cell  Cultures 

The  human  thyroid  epithelial  cell  line  (HTori-3) 
used  in  our  study  was  originally  derived  by  trans¬ 
fecting  the  epithelial  cells  with  a  plasmid  contain¬ 
ing  the  origin  of  replication  defective  SV40  ge¬ 
nome  [14],  The  cells  were  cultured  in  medium 
consisting  of  equal  volumes  of  Dulbecco's  modi¬ 
fied  Eagle’s  medium  and  Ham's  FI 2  medium 
supplemented  with  lck  fetal  calf  serum,  penicillin 
(10  units  per  ml),  streptomycin  ( 100  pg  per  ml)  and 
L-glutamine  (2  mM).  It  has  previously  been  dem¬ 
onstrated  that  following  irradiation  of  these  cells  in 
vitro  with  either  single  doses  of  gamma  irradiation 
or  alpha  particle  irradiation  that  tumors  developed 
in  athymic  nude  mice  following  transplantation  of 
the  irradiated  cells  [4,6.7],  Cells  were  not  pooled  at 
any  stage  so  that  the  resulting  tumors  must  be  re¬ 
garded  as  independent  events.  The  unirradiated 
cells  in  these  series  of  experiments  did  not  produce 
tumors  in  athymic  nude  mice.  Cell  lines  were  de¬ 
rived  from  the  primary  tumors.  Three  of  these  tu¬ 
mor  cell  lines  will  be  utilised  in  this  study.  The  first 
line  (HTG1)  was  derived  from  cells  exposed  to 
0.5  Gy  gamma  irradiation.  A  tumor  was  detected 
17  weeks  after  injection  into  athymic  nude  mice 
and  was  removed  at  25  weeks  when  the  tumor  was 
8  mm  in  diameter.  A  further  two  tumor  cell  lines 
(HTA3  and  HTA6)  were  derived  from  cells  ex¬ 
posed  to  0.25  Gy  and  0.5  Gy  alpha  particle  irradia¬ 
tion.  The  tumor  HTA3  was  detected  at  13  weeks 
after  injection  and  was  removed  at  1 7  weeks  when 
the  tumor  diameter  was  8  mm.  Similarly  HTA6  was 
detected  at  17  weeks  and  removed  at  23  weeks 
when  two  tumors  of  6  mm  in  diameter  were  pres- 
nt.  Following  establishment  of  cell  lines  in  vitro 
from  the  primary  tumors,  cloned  cell  lines  were 
derived.  From  the  HTG1  primary  cell  line,  two 
lines  were  derived  (HTGlcll  and  HTGlcl3).  Simi¬ 
larly  cloned  cell  lines  were  derived  from  the  HTA3 
primary  cell  line  (HTA3cl3)  and  the  HTA6  primary 
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Fig.  1.  HTori  cell  lines  established  from  primars  thyroid 
epithelial  cells  following  SV40  immortalization  and  radia¬ 
tion-induced  transformation.  Cell  lines  had  been  tested  lor 
tumorigenicity  in  athymie  nude  mice.  Cultures  were  estab¬ 
lished  from  the  primary  tumors  that  developed  and  cell  lines 
generated  from  cloned  cells. 

cell  line  (HTA6cl2,  HTA6cl4,  HTA6c)5.  HTA6cl6, 
HTA6cl7  and  HTA6cl9).  The  origins  of  the  cell 
lines  are  summarised  in  Figure  1. 

Tumorigenicity  Screening 

The  cloned  cell  lines  derived  from  primary  tumor 
cell  lines  HTG1  (two  lines),  HTA3  (one  line)  and 
HTA6  (six  lines)  were  screened  for  tumorigenicity 
in  athymie  nude  mice.  Cells  (2-3  x  106)  were  in¬ 
jected  subcutaneously  into  the  dorsal  region  of 
three  athymie  nude  mice.  The  mice  were  screened 
for  tumor  development  for  5  months  and  scored 
positively  if  tumors  of  greater  than  3  mm  diameter 
growing  progressively  were  detected. 

G-Banding 

Cells  were  seeded  onto  glass  slides  and  metaphase 
spreads  were  prepared  after  in  vitro  culture  of  cells 
for  several  days  following  standard  procedures.  Af¬ 
ter  G-banding  with  Wright’s  staining  solution. 


Tumor  Cells 

karyotyping  was  performed  according  to  the  Inter¬ 
national  System  for  Human  Cytogenetic  Nomen¬ 
clature  (1SCN)  1 15]. 

Spectral  Karyotyping 

SKY  analysis  was  performed  as  described  previ¬ 
ously  [10].  Briefly,  chromosome  preparations 
were  aged  for  two  weeks  at  room  temperature. 
Metaphase  spreads  were  pretreated  with  RNase 
A  (0.1  mg/ml  in  2  x  SSC)  and  pepsin  solution 
(10  pg/ml  in  0.01  M  HC1)  and  fixed  in  1%  form¬ 
aldehyde.  Pepsin  digestion  was  performed  under 
microscopic  control.  Slides  were  placed  in  denatur¬ 
ing  solution  (70%  formamide  | FA]  in  2  x  SSC)  at 
70°C  for  150  sec  and  dehydrated  in  a  70%,  80'/, . 
100%  ethanol  series.  Subsequently,  spreads  were 
hybridized  with  a  probe  mixture  supplied  by  Ap¬ 
plied  Spectral  Imaging  (AS1.  Carlsbad,  Calil.).  The 
probe  mixture  contained  24  whole  chromosome 
painting  probes,  labeled  with  combinations  of  five 
different  fluorescent  dyes  (Spectrum  Green.  Spec¬ 
trum  Orange,  Texas  Red,  Cyanine  ICy]  5  and 
Cy5.5)  [16].  The  probe  (SKY™  mixture)  was 
denatured  at  75°C  for  7  min  and  incubated  at  37°C 
for  1  hr.  The  hybridization  solution  was  applied  to 
the  denatured  metaphases  and  incubated  lor  two 
days  at  37°C.  Posthybridization  washes  were  per¬ 
formed  in  50%  FA/2  x  SSC  (three  times  for  5  min) 
at  45°C,  1  x  SSC  (twice  for  5  min)  and  in  4  x 
SSC/0.1%  Tween  20,  (2  min  at  20°C).  Detection 
of  biotinylated  probes  was  performed  with  avidin- 
Cy5.  Digoxigenin-labeled  probes  were  detected 
with  a  mouse  anti-digoxigenin  antibody  followed 
by  a  goat  anti-mouse  antibody  conjugated  to 
Cy5.5.  Metaphase  spreads  were  stained  with 
4'.6-diamidino-2-phenylindole  (DAP1)  solution 
(150  ng/ml  in  2  x  SSC)  and  covered  with  antifade 
solution  (Vectashield  mounting  medium;  Vector 
Laboratories,  Burlingame,  Calif.).  Metaphase  im¬ 
ages  were  acquired  by  use  of  a  SpectraCube  system 
(AS!)  and  analyzed  with  the  SKYView  imaging 
software  [17].  Aberrations  were  described  accord¬ 
ing  to  the  1SCN  nomenclature  [15]. 

Comparative  Genomic 
Hybridization  (CGH) 

CGH  was  performed  on  cell  lines  as  described 
[18,19]  with  slight  modifications.  Briefly,  whole 
genomic  DNA  was  isolated  from  the  cell  lines  ac¬ 
cording  to  standard  procedures  and  labeled  with 
biotin- 16-dUTP  (Roche  Molecular  Biochemicals, 
Mannheim,  Germany)  and  biotin- 1 1-dCTP  (Life 
Technologies,  Eggenstein,  Germany)  by  nick- 
translation.  Normal  female  reference  DNA  was 
isolated  from  peripheral  lymphocytes  of  a  healthy 
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donor  and  labeled  with  cy3-eonjugated  dUTP  and 
dCTP  (Amersham.  Braunschweig.  Germany).  After 
hybridization  to  normal  metaphase  spreads  of  a 
healthy  male  donor,  the  biotin-labeled  DNA  was 
detected  with  avidin-FITC  (Vector  Laboratories) 
(diluted  in  0.1  M  sodium  phosphate,  pH  8.0,  0.1% 
nonidet-P40.  0.1%  sodium  azide  plus  5%  nonfat 
dry  milk).  After  washes  in  PN  buffer  (0.1  M  so¬ 
dium  phosphate.  pH  8.0.  0.1%  nonidet-P40),  slides 
were  counterstained  with  DAP1  (150  ng/ml  in 
2xSSC)  for  karyotyping  and  mounted  in  antifade 
solution.  For  CGH  analysis,  at  least  10  metaphases 
were  imaged  and  karyotyped  using  a  Zeiss  Ax- 
ioplan  2  fluorescence  microscope  equipped  with 
filter  sets  (single-band  excitation  filters)  for  DAPI, 
fluorescein,  and  tetramethylrhodamine  isothiocya¬ 
nate  (TR1TC ).  Averaged  profiles  were  generated  by 
a  CGH  analysis  software  (MetaSystems,  Al- 
tlussheim.  Germany)  from  at  least  10  to  15  homolo¬ 
gous  chromosomes  and  interpreted  according  to 
published  criteria  1 18.20]. 

RESULTS 

Ten  HTori  cell  lines  were  investigated  cytogeneti¬ 
cally  by  G-banding.  SKY.  and  CGH,  among  them 
the  parental  cell  line  HTori-3  as  well  as  nine  radia¬ 
tion-induced  tumor  cell  lines  (Fig.  1 ).  The  results  of 
the  tumorigenieity  screening  of  the  nine  cloned  cell 
lines  are  summarised  in  Figure  1.  Most  of  the 
cloned  lines  were  tumorigenic  in  athymic  nude 
mice.  Two  of  the  cloned  lines  (HTA6cl4  and 
HTA6cl9  )  derived  from  the  HTA6  primary  cell 
line  failed  to  produce  tumors. 

G-banding  analysis  revealed  a  near-pentaploid 
chromosome  number  in  cell  line  HTA3cI3,  while 
all  other  cell  lines  appeared  to  have  a  hypotriploid 
chromosome  number  and  various  marker  chromo¬ 
somes  that  could  not  be  classified  further. 

Therefore,  a  detailed  cytogenetic  characteriza¬ 
tion  was  only  feasible  combining  results  from  SKY 
and  CGH  analysis.  For  SKY  analysis,  about  4  to  10 
metaphases  were  analyzed  from  each  cell  line.  A 
total  of  79  clonal  markers  (Ml  to  M79),  i.e.,  struc¬ 
turally  rearranged  or  deleted  chromosomes  appear¬ 
ing  in  at  least  two  metaphases,  were  identified 
(Table  1).  Some  of  the  markers  have  very  similar 
breakpoints  (M1/M2,  M9/M11,  M20/M21,  M23/ 
M24,  M44/M47.  M58/M59/M60,  M69/M70)  and 
differ  only  in  one  translocation  partner  chromo¬ 
some. 

Twenty-three  different  immortalization  mark¬ 
ers  were  identified  by  SKY  analysis  of  the  parental 
cell  line.  They  are  indicated  in  Tables  1  and  2  with 
a  single  underline  and  examples  of  some  of  these 
markers  are  shown  in  Figure  2.  Additional  changes 


in  the  immortalized  HTori-3  cell  line  arc  gains  on 
4p.  5q,  6p.  16p.  1 7q.  20,  21,  and  22  as  well  as 
losses  on  3p.  1  Iq.  I8q.  and  X  (Table  3).  Cytoge¬ 
netic  markers  that  appeared  exclusively  in  cell  lines 
originating  from  y-my-induced  tumors  (HTGlcll, 
HTGlcl3)  or  in  cell  lines  developed  from  a-ray- 
induced  tumors  (HTA3cl3  and  HTA6  clones)  are 
marked  in  Table  2  by  dotted  and  double  underlines, 
respectively.  In  cell  lines  originating  from  y-ray- 
induced  tumors,  markers  Ml 6,  M26,  M27,  M30, 
M35.  M39.  M41.  and  M52  became  apparent, 
whereas  markers  M6.  M8.  MI 5,  Ml 7,  Ml 8,  M21, 
M34.  M69.  and  M72  were  detected  exclusively  in 
cell  lines  developed  from  a-ray-induced  tumors 
(Fig.  3). 

Cytogenetic  markers  obviously  related  to  the 
tumorigenic  process  were  studied  in  HTA6  clones 
showing  a  different  tumorigenic  behavior.  Clones 
HTA6cl2  and  HTA6cl6  induced  the  same  tumor 
frequency  as  the  corresponding  primary  culture, 
whereas  HTA6cl5  and  HTA6cl7  showed  a  reduced 
and  HTA6cl4  and  HTA6cl9  a  completely  missing 
tumorigenic  behavior  in  nude  mice.  In  the  tumori¬ 
genic  lines  HTA6cl2.  HTA6cl6,  HTGlcll.  and 
HTGlcl3.  loss  ot  9q32-q34  became  apparent.  Ad¬ 
ditionally.  loss  of  7q2l-q3i  and  an  increased  copy 
number  of  chromosome  20  (four  copies)  appeared 
to  be  indicative  for  tumorigenic  lines.  The  nontu- 
morigenic  cell  line  HTA6cl9  as  well  as  lines 
HTA6cl5  and  HTA6cl7  with  a  reduced  tumorigenic 
potential  exhibited  only  two  copies  of  chromosome 
20  (Tables  2  and  3). 

Besides  these  differences,  common  chromo¬ 
somal  changes  were  observed  in  cell  lines  originat¬ 
ing  from  different  irradiation  experiments.  Loss  of 
13cen-q21  occurred  in  each  irradiated  cell  line  and, 
thus,  independently  in  three  different  irradiation  ex¬ 
periments.  Loss  of  8q23  as  well  as  gain  of  lq32- 
qter  and  2q  1 1  2-q  14. 1  were  visible  in  HTG1  and 
HTA6  clones  (Table  3). 

DISCUSSION 

The  malignant  transformation  of  human  epithelial 
cells  in  vitro  by  the  combined  action  of  the  DNA 
tumor  virus  SV40  and  ionizing  y-  and  a-radiation 
is  paralleled  by  cytogenetic  changes.  Previous  re¬ 
ports  of  in  vitro  systems  developed  to  study  single 
steps  of  radiation-induced  neoplastic  transforma¬ 
tion  of  epithelial  cells  1 1.2.5,21]  already  postulated 
a  multistep  process  for  cancer  development.  It  is 
associated  with  chromosomal  alterations  that  were 
also  detected  in  our  in  vitro  system  to  investigate 
the  radiation-induced  transformation  of  the  thyroid 
epithelium  [4,6.7,14],  To  search  for  specific 
changes  attributable  to  immortalization,  radiation 
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Table  1.  Clonal  Marker  Chromosomes  in  HTori  Cell  Lines  Analyzed  by  SKY 


Marker  Chromosomal  composition _ Marker 


Chromosomal  composition 


Ml  der(  I )( I pter  »q44::5q3I  -qter) 

M2  dcr(  1 )( 1  pier  >q44::9ql  I  -qter) 

M3  der(  1  )(20pter  *p  1 1 : :  I p V 1  -qter) 

M4  der(  I )( 1  pier  -q 1 1  .various  chromosomes) 

M3  der(  1 )( Iqter  ►pi  l::various  chromosomes) 

Mh  der(  I )( 1  q2 1  ►pier::)5qter  -q22) 

M7  del(2)(:p23  -qter) 

MX  i(2)(qter  ♦ql()::qlO  -qter) 

M9  der(2)(2qter  -pll::3ql  I  -q26::7p?q?> 

M 10  der(2)(2pter  -q  1 1  ::6q  16  -qter) 

MU  der(2)(2qter  -pl  l::10p?q?) 

M12  der(3)(7qier  -q22::3p25  -q25::4pl3::3q25  -qter) 

M13  der(3)(8qter  >q?::3pll  -qter) 

M14  der(3)(  Iqter  -q32::l2q!4  -q  1 2::  1 9p?q?: :2q 1 4  ► 

.  qll.2::3pll  -qter) 

M 13  der(4)(4p?q?::4p!6  -q32:) 

M 1 6  der(4)(  I  Oqter  -q  1 1 :  :4p  1 5  -qter) 

M 1 7  der(4){  20qter  ►q  1 1 :  :4p  1 4  -qter) 

MIX  del(5)(pter  -ql  1:) 

M 19  der(5)(3qter  -q21  ::5pl5  -qter) 

M20  der(5)(5qter-pl5::8q?::14pl  I  ~-qter) 

M2l  der(5)(5qter  ►pl5::8q?::6q)5  -q23::8q?::14pl  1  -qter) 

M22  del(6)(pter  -ql3:) 

M2 3  del(6)(pter -q15:) 

M24  der(6)(pter-q  1 5: :  1 2p?q?) 

M25  der(6)(6pter  -q21::2p21  "pter) 

M26  der(6)(9pter  -pll::6q1I  -pter) 

M27  der(6)(6pter-q27::16ql  1  -qter) 

M28  der(7)(7pter  -q32 : :  I Oq 1 1  -qter) 

M29  der(7)(7pter-ql  1  ::20p  1 2  -pter) 

M30  der(7)(7pier  -ql  1  ::22ql  I  -qter) 

M3 1  i{8)(qter-ql0::ql0  -qter) 

M32  der{9)(9pter-q33::  1 3q22  -qter) 

M33  der(9)(21qter  -ql  1  ::9pl  1  -qter) 

M34  der(9)(9qter  -pl  1 : :  1 5q  1 1  -qter) 

M35  der(10)(4pter-pl5::10pll  -pter) 

M36  del(  1 1  )(pter~-q2 1 : ) 

M37  der(  1 1 )(:  1  Ipl  1  --qter) 

M38  i(  1 1  )(qtei  -q  1 0: :q 1 0 ^qter) 

M39  died  1  ;1 9)(  1 9qter  -p  1 3 : :  11  p  1 1  --qter ) 

M40  der(  11 )( 1 1  pter->q  1 3::Xq  1  l~-qter) _ 


M4I  den  I2)i  I2pter  -q!3::20p1 1  -pier) 

M42  der(  1 3 )( 9pter  -pll::l3pl  1  --qter) 

M43  diet  14:20)1 2()pter  ►qll::I4p?  -qter) 

M44  derl  14 )( 2 Iqter  ►qll::14pl  1-qter) 

M43  deri  14)(  I4qter  ►p?:: Ipl  1-pter) 

M46  dert  14 )( Iqter  -lql  1 ::  1 4p?  ►!4q32::17q21  ►qter) 

M47  den  1 3)t2 Iqter  -ql  l::15pl  I  -qter) 

M48  der(  13)(  I7pter  -pi  l::15p?  -qter) 

M49  dert  1 5 )( 7pter  ►pll  ::15p?  -qter) 

M50  del(16)(pter  ►qll:) 

MSI  denl6)(16qter  ►16pl3::6p22  -pter) 

M52  den  I7)(  I7pter  -q23::3p‘?q?) 

M53  den  1 7 )( 17qter  -pl  l::8q?::l  Iql4  -q23::Xp?q?:: 1 9p?q?::  1 6p 1 2  ►pter) 
M54  dert  1 7  K  1 1  qter  -q!3::17pl3  *q  1 1 :  :8q?  1 1  q  1 4  ► 
q23::Xp?q?::  1 9p?q?::  1 6pl  2  ►pter) 

M33  it  1 8 Xqier  -qlO::qlO  -qter) 

M56  delt  18)( pter  -p21:) 

M57  delt  1 8){  :pl  1  -qter) 

M58  derl  I9)(  I9pter  -q!3::lq?) 

M59  dert  19)(  19pter  ►q!3::5p?q?) 

M60  dert  19)(  19pter  -q!3::l  lql 3 --qter) 

M61  dert  19)(  19qter  ►p 1 3::22q  1 1  -q 1 2::  1 7q22  -qter) 

M62  der(20)(20pter  -ql  1 : :  1 2q  1 1  --12qter) 

M63  dert  20>(  19pter  ►qll  ::20ql  I-pter) 

M64  der(20)(20pter  -q  1 1 ::3p?q?::2 1  q 1 1  -qter) 

M63  der(20)l20qter  -p  1 3: :  1 6p?q?::  1 4q  1 1  -qter) 

M66  der(20)(:17q23  -ql  1  ::20ql  1  -p!3:: 16p?q?:: I4q  -qter) 

Mn7  dert  20 )( 2()qter  *20pl  3::5p?q?::14ql  I  *qter) 

M68  it 21  Xqier  ►q IO::qlO  -qter) 

M69  dert 2 1 )( 2 Iqter  ►pl  1  ::5ql2-5qter) 

M70  dert 2 1 )( 2 Iqter  ►pl  1 ::  !0q22-qter) 

M71  dert 21  )(20pter  -pl  1  ::2 Ipl  1  “-qter) 

M72  derl 21 )( 17 pter  ►pi  l::21pl  1-qter) 

M73  der(2l :22(2 Iqter  -pi  1  ::22pl  1 --q  1 2::  1 7q22  -pter) 

M  74  dert  22 )( 22qter  -p  1 1 : :  1  Op?q?: :  1 8q  1 1  -qter) 

M75  dert 22 ){ 22p 1 3  -qI2::17q22-qter) 

M76  der(22)(  I7qter  -q22::22ql2^pl  l::9q!2  -q33::  13q22  -qter) 

M77  der(22)(22qter  ►pll ::18ql  1-qter) 

M78  der(22)(  17qter  -q22::22ql2~-pl  1::1  lql  I  -qter) 

M79  dert XX  lOpter  -pH  ::Xql  1-pter) 


Bold  and  single  underline,  markers  of  SV40  immortalization. 


effects,  and/or  tumorigenicity,  we  compared  chro- 
mosomal  markers  among  the  cell  lines  presented  in 
Figure  1.  The  parental  cell  line  HTori-3  [14]  with 
cytogenetic  alterations  obviously  related  to  SV40 
immortalization  served  as  starting  point  for  a  com¬ 
parison  with  radiation-induced  tumorigenic  cell 
lines.  Primary  cultures  from  tumors  developed  in 
athymic  nude  mice  (HTG1,  HTA3,  HTA6;  Fig.  1) 
provided  information  on  chromosomal  changes  in¬ 
duced  by  different  doses  of  two  radiation  qualities 
in  addition  to  alterations  that  developed  during  the 
process  leading  to  tumorigenicity.  To  identify  ab¬ 
errations  that  are  indispensable  for  the  tumorigenic 
process,  clones  from  the  primary  cultures  of  the 
tumors  were  rechecked  for  tumorigenicity 


(HTGlcll,  HTG lcl3,  HTA3cl3,  HTA6cl2, 
HTA6cl4,  HTA6cl5,  HTA6cl6,  HTA6cl7, 
HTA6cl9;  Fig.  1). 

The  G-banding  analysis  of  the  parental  cell 
line  HTori-3  suggested  that  polyploidization  and 
marker  chromosomes  already  originated  during  the 
SV40-mediated  immortalization  step.  Karyotypic 
instability  after  SV40  integration  preceding  neo¬ 
plastic  transformation  has  previously  been  demon¬ 
strated  for  human  fibroblasts  [22],  Recurrent  and 
more  likely  important  markers  for  immortalization 
are  represented  by  Ml,  M2,  M12,  M14,  M20,  M32, 
M36,  M40,  M51,  M53,  and  M79  because  they  oc¬ 
curred  in  almost  every  cell  line  investigated  (Table 
2,  Fig.  2).  In  addition,  CGH  revealed  numerous 
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Table  2.  Clonal  Karyotypes  of  HTori  Cell  Lines  Analyzed  by  SKY 


HTori-3 

HTGl  -clone) 

HT(i  l-clone3 

HTA3-clone3 

Chromosome  number 

67-70 

62-65 

62-67 

97-122 

Ploitly  ranee 

3n-’ 

3n-;i 

3n-“ 

5n+h 

Chromosome  1 

1  M1;M2 

1  :M1;M2 

1 : 1  :M1:M2 

1:!:M1;M6 

Chromosome  2 

2:2;2:M10;\11J 

2:2:2 

">.  s .  s 

2:2:M8;M9 

Chromosome  3 

3;3;M12;M14 

3:M12;M14 

3:\112;M  14 

3:3  M12;M13 

Chromosome  4 

4:4:4 

4:4:M  16 

4:4:M  16 

4:4 :4:4:4 

Chromosome  5 

5:5:M20 

5:5:M20 

5:5:M20 

5:5  :M  1 8:  M 1 9;M20 

Chromosome  6 

6:6: M 24 

6:6:M26;M27 

6:6:M26:M27 

6:6:6:M22;M25 

Chromosome  7 

7;7:7:M28 

7:7;M30 

7:7:M30 

7:7;7:7:7:7 

Chromosome  X 

X:8 

8:8 

X;X 

X:8:8:8:8:M31 

Chromosome  9 

9;9:M32 

9:M32 

9:M32 

9:9:M32;M34 

Chromosome  10 

10: 10;  10 

10:1 0:\135 

10:10:M35 

10:10:10 

Chromosome  1  1 

1  1  :M36;M37;M40 

M36:M39;M40 

M36;M39:M40 

11:11  :M36;M38 

Chromosome  12 

12:12:12 

12: 12:M41 

1 2;  1 2:M41 

12:12:12:12 

Chromosome  13 

13:13:13 

13:13 

13:13 

13:13:13:1 3:M42 

Chromosome  14 

14:14 

14:14 

14:14 

14:14 

Chromosome  15 

15:15:15 

15:15:M47 

15:15:15 

15:15 

Chromosome  16 

16:M50;M51 

I6:M50;M51 

I6:M50:M51 

16;  16:16:M51 

Chromosome  17 

17:M53 

17:M52M53 

17:M52:M53 

17:17:17  :M53 

Chromosome  IS 

IX:  18:M56 

18:18:M56 

1X:M56 

18:18:18 

Chromosome  19 

I9:M58;M59 

19 

19 

19:19:19:M60 

Chromosome  20 

20:20:20:20 

20:20:20:20:20;M64 

20:20:20:20:20 

20;20;20:20:20:20;M62 

Chromosome  21 

21:21  ;2!  :M68,M70 

21:21  :M68 

21:21  :M7 1 

21:21:21  :M69 

Chromosome  22 

22:22:22:M75 

22:22:M75 

22:22  :M75 

22:22:22 

X  chromosome 

X:X:M79 

X:M79 

X:M79 

X:X 

chromosomal  imbalances  in  HTori-3  (Table  3). 
which  might  be  attributed  to  the  immortalization 
process.  Of  special  interest  is  the  finding  of  loss  of 
6q  material.  SEN6.  a  locus  for  SV40-mediated  im¬ 
mortalization  of  human  cells,  maps  to  6q26-q27. 
and  has  been  demonstrated  to  be  deleted  in  SV40 
immortalized  cells  [23].  Experiments  using  micro¬ 
cell-mediated  chromosome  transfer  demonstrated 
that  chromosomes  1.2.4.  6.  7.  and  1 1  can  confer 
senescence  in  various  tumor  cell  lines  [24-30].  It  is 
further  known  that  chromosome  llq  is  a  strong 
candidate  region  for  proto-oncogenes  playing  a  key 
role  in  immortalization  of  keratinocytes  upon  trans¬ 
fection  with  human  papillomavirus  [31].  In  our  in 
vitro  model  system,  chromosomes  lq,  4,  6,  and  1 1 
were  found  frequently  affected  by  translocations  as 
well  as  by  numerical  imbalances  (Fig.  2,  Table  3), 
suggesting  that  such  changes  occur  during  immor¬ 
talization  preceding  tumorigenicity.  It  was  assumed 
that  inactivation  due  to  chromosomal  alterations  of 
genes  controlling  cellular  senescence  can  activate 
the  process  of  immortalization  [32].  Our  SKY  and 
CGH  data  on  HTori-3  cells  demonstrated  that  other 
immortalization-related,  so  far  nonreported  aberra¬ 
tions  affected  chromosomes  3p,  8p,  9p,  16p,  19q, 
21,  22,  which  might  harbor  genes  important  for  the 
immortalization  process.  Thus,  this  data  can  serve 
as  a  starting  point  for  further  microcell-mediated 


assays  to  investigate  their  role  in  senescence  and 
immortalization. 

The  original  tumor  cell  lines  (HTG1.  HTA3, 
HTA6)  exhibited  extensive  chromosomal  instabil¬ 
ity  (data  not  shown),  which  was  less  prominent  in 
clones  derived  from  these  cell  lines.  We  therefore 
decided  to  interpret  only  data  originating  from  the 
cloned  cell  lines.  Whereas  markers  Ml 6,  M26, 
M27,  M30.  M35.  M39.  M41,  and  M52  became 
apparent  only  in  the  y-ray-induced  lines  (HTGlcll 
and  HTG)cl3).  markers  M6.  M8.  Ml 5,  Ml 7,  Ml 8, 
M21,  M34,  M69,  and  M72  were  detected  only  in 
the  a-ray-induced  lines  (HTA3cl3  and  HTA6 
clones).  A  similar  observation  was  made  by  CGH 
analysis,  which  revealed  losses  on  lq  and  2  as  well 
as  gains  on  3q  and  9p  in  HTA3cl3  only.  If 
HTA3cl3  and  HTA6  clones  were  compared  follow¬ 
ing  SKY  analysis,  it  was  observed  that  M6,  M8, 
Ml 8,  M34,  and  M69  appear  only  in  HTA3cl3, 
while  Ml 5,  Ml 7,  M21,  and  M72  occurred  exclu¬ 
sively  in  HTA6  clones.  The  copy  numbers  for  chro¬ 
mosomes  7,  8,  and  12  were  also  different  between 
HTA3cI3  and  HTA6  clones.  If  clones  1  and  3  of 
HTG1  cell  line  were  compared,  only  small  differ¬ 
ences  became  apparent,  probably  reflecting  a  clonal 
evolution  of  tumor  cells.  Also,  HTA6  clones  2, 4,  5, 
6.  7,  and  9  were  very  similar  in  their  aberration 
pattern,  although  there  were  more  extensive  alter- 


172  Zitzelshergcr:  Radiation-Transformed  Thyroul  Tumor  Cells 
Table  2.  Continued  _ 


HTA(vdone2  HTA6-cloned4  HTA6-clone5 


Chromosome 
Ploidy  range 
Chromosome  I 
Chromosome  2 
Chromosome  3 
Chromosome  4 
Chromosome  3 
Chromosome  6 
Chromosome  7 
Chromosome  8 
Chromosome  9 
Chromosome  10 
Chromosome  1 1 
Chromosome  12 
Chromosome  1 3 
Chromosome  14 
Chromosome  15 
Chromosome  10 
Chromosome  17 
Chromosome  IX 
Chromosome  1 9 
Chromosome  20 
Chromosome  2\ 
Chromosome  22 
X  chromosome 


59-64 
3n-  1 
l.M  l:\12 


M12;M  14 


59-  6X 
3n-:i 

M2;M3:M4;M5 

3:M12;M14 


57-64 

3n-u 

1;M1;M2 

2:2:2 

3:M12;M14 


,M15;MI7  4;M15;M17  4;M15;IVn7 


.5:M21 

M22:M23 

:7 

,8 

>132 

10:10 

M36;M37;M40 

12:12 

13:13 


14:14 

15:15 

16:16 


>144 

15 

M50:M51 


I7:M53 
1  X;M56 
19 


5:5:M21 

6:6;M23 

7:7;7:7 

X:X 

9:M32;M33 

10:10:10 

M36:M40 

12:12 

13:13 

M45:M46 

15:M48;M49 

16:M50;M51 

1 7:M54 

1X:1X:M56 

19 


20:20: 

21:21 

“f 


:20:20:M63 

>168;M72 

;M75 


X:M79 


20;20:M63 
M72 
M75 
X:M79 


20:20 

21:21 


HTA6-clone6 


HTA6-clone7  HTA6-clone9 


5:5:M21 
6:6:M23 
7:7:M29 
8:8 

9:M32 
10:10:10 
M36;M37;M40 

12:12 

13:13 

14 

15:15:15 

16:M50,M51 
1 7:M53 
18 

19:M61 

20;20;M63;M67 

21:21:M68;M72 

22:M75;M76;M77 

X:M79 


59-64 

3n-‘l 


1  :M1;M2 


3:M12;M14 
4:M15;M17 
5:5:M18:M21 
6:6:M23 
7:7:M29 
8:8 

9:M32 
10:10 

M36;M37;M40 

12:12 

13:13 

14:14 

15:15:15 

16:M50;M51 

17:M53 

18:M56;!VI57 

19 

20:20:20:2():M63;M67 

21:21:M68;M72 

22:22  :M73;M75;M76;M  78 

X:M79 


58-61 

3n-“ 

1  :MI  ;M2 

2:2:M7 

3:M12;M14 


60-63 

3n-il 

I:M1;M2 

2:2:2 

3:M12;M14 


4:M  I5;M17  4:M15:M17 


5:5:M21 
6:6:M22 
7:7:M29 
8:8 

9:M32 

10:10:10 

M36;M37;M40 

12:12 

13:13 

14:14 

15:15:15 

16:M50;M51 

17:M53 

18 

19 

20:20:M63;M66 

21:21:M68;M72 

22:M74:M76 

X:M79 


5:5:M21 
6:6:1V!  23 
7:7:M29 
8:8 

9.M32 
10:10:10 
M36;M37;M40 
12:12 
13:13 
!4;M43 
15:15:15 
16:M50;M51 
1 7:M53 
18:M55;M56 
19 

20;20:M63;M65 
2 1  :M68:M72 
22:M76 
X:M79 


Bold,  markei  ehromoMimes.  single  underline,  markers  ol  SV40  immortalization:  dotted  underline,  markers  characieristn 
markers  ehaiaeteiistie  toi  o-iav-indueed  lines 
Mil-,  hypotriplonl  veil  line 
h5n±.  near-peniaploul  eell  line 


lor  'v-riiN  -  induced  line:  double  underline. 


ations  than  between  the  HTG  clones.  Some  of  these 
different  aberrations  might  be  responsible  for  the 
varying  tumori genic  behavior  of  the  HTA6  clones. 
Moreover,  it  appeared  that  deletions  on  chromo¬ 
somes  9q32-q34  and  7q21-q31  as  well  as  increased 
copy  numbers  of  chromosome  20  were  important 
for  the  tumorigenic  phenotype.  According  to  our 
data,  it  seems  that  varying  combinations  of  these 
three  aberrations  could  contribute  to  a  more  aggres¬ 
sive  tumorigenic  behavior.  This  is  concordant  with 
findings  in  a  variety  of  solid-tumor  entities  that 
frequently  exhibit  losses  on  7q21-q31  and  9q32- 
q34  as  well  as  a  gain  of  chromosome  20  [33-35].  It 
can  be  further  assumed  that  common  aberrations 
among  radiation-induced  lines  (loss  on  8q23  and 
13cen-q21 .  gain  on  lq32-qter  and  2ql  1.2-ql4.1) 
are  strong  indicators  for  the  involvement  of  genes 
located  within  these  regions  in  radiation-induced 
tumor  development  in  thyroid  epithelial  cells.  Al¬ 
though  they  probably  imply  significant  genetic 
changes,  it  is  likely  that  they  were  not  sufficient  to 
cause  the  aggressive  malignant  phenotype  since 
they  occurred  also  in  the  nontumorigenic  lines.  It 
is,  however,  of  interest  that  recent  CGH  and  loss  of 
heterozygosity  (LOH)  studies  on  follicular  carcino¬ 
mas  |36,37]  and  Hurthle  cell  tumors  [38]  of  the 


thyroid  gland  revealed  the  same  aberrations  (losses 
on  chromosomes  1 3q  and  9q3 1  -q33.  gains  on  chro¬ 
mosomes  lq  and  20)  as  this  model  cell  system  of 
radiation-induced  transformation  of  epithelial  thy¬ 
roid  cells. 

A  comprehensive  breakpoint  analysis  of  chro¬ 
mosomal  aberrations  in  HTori  lines  (Fig.  4)  and  in 
radiation-induced  childhood  thyroid  tumors  from 
Belarus  [10]  revealed  coinciding  mapped  break¬ 
points  on  chromosomes  lq,  4q,  5q,  6,  7q,  9q,  lOq, 
1  lq,  12q,  15q,  16,  and  22q.  However,  some  clus¬ 
tered  chromosomal  breakpoints  observed  in  HTori 
lines  on  chromosome  21  and  most  prevalently  on 
chromosome  20  were  different  from  the  Belarus¬ 
sian  tumors.  These  findings  indicate  that  radiation- 
induced  thyroid  tumors  might  share  some  hot  spots 
of  chromosomal  damage.  Different  doses  and  ra¬ 
diation  qualities  may  have  led  to  variations  in  the 
aberration  pattern. 

In  conclusion,  we  demonstrate  specific  chro¬ 
mosomal  aberrations  associated  with  SV40  immor¬ 
talization  and  radiation-induced  transformation  of 
human  epithelial  thyroid  cells.  The  data  are  in  line 
with  previous  reports  describing  immortalization- 
related  changes  in  epithelial  cells  but  also  indicate 
newly  identified  chromosomal  markers  originating 


M20 


M32 


M36 


M40 


M51  M53  M79 


Fig.  2.  Chromosomal  markers  of  SV40  immortalization  detected  by  SKY.  Markers  Ml.  M2.  MI2,  M14,  M20,  M32,  M36. 
M40.  M3 1 .  M53.  and  M79  are  shown  as  RGB  image.  DAPI-banded  image,  and  classified  image  in  pseudocolors  (from  left 
to  right)  The  structural  composition  of  each  marker  chromosome  is  indicated  in  Table  I. 

A  HTA3  cell  lines 


M6  M8  M18  M34  M69 


B  HTA6  cell  lines 


M15  M17  M21  M72 


C  HTG1  cell  lines 


M16  M26  M27  M30  M35  M39  M41  M52 

Fig.  3.  Chromosomal  markers  of  different  HTori  cell  lines.  Specific  markers  of  HTA3  clones  (A).  HTA6  clones  (B),  and 
HTGI  clones  (C)  are  demonstrated  as  RGB  image.  DAPI-banded  image,  and  classified  image  in  pseudocolors  (from  left  to 
right).  The  chromosomal  composition  of  each  marker  is  indicated  in  Table  I. 


Fig.  4.  Breakpoints  of  structural  chromosome  aberrations  and  chromosomal  gains  and  losses  in  different  HTori  cell  lines. 
They  were  indicated  for  the  individual  cell  lines  in  different  colors  (blue,  HTon-3;  green,  HTG1;  grey,  HTA3;  red,  HTA6). 
Breakpoints  of  structural  aberrations  are  shown  on  the  left  side  of  each  ideogram.  On  the  right  side  of  each  ideogram,  gains 
are  represented  by  running  lines  and  losses  by  dashed  lines. 
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Tabic  3.  CGH  Results  of  10  HTori  Cell  Lines 

Cell  line  Gains  Losses 


HTori- 3 


HI  G I - 
clone  I 
I II G  I  - 

clone 3 
HTA3* 
elone3 


HTAft 
cl  one  2 

HTAft- 

clone4 

HTAft- 

clom*5 

HTAft- 

cloneft 

HTAft- 

clone7 

HT  Aft- 
clone*) 


4p15,  5,  6p12-pter,  7,  8q24.1-qter,  9q21-q22,  10,  1  lq!3, 
I2pter-ql5,  13q22-qter,  14q24-qter,  15q22-qter, 
lftpl2-pter.  )7q21-qter,  19p,  20  (a:  20q13>,  21,  22 
lq32-qter.  2ql  I.2-ql4.1.  4pl5,  5q31-qter,  6p,  lOp, 
12ql3-ql5,  16pl2-pler,  17q2Lqter,  20  7a:  20pl3).  21,  22 
Iq32-qter.  2ql  1  2-ql4.  4pl5,  5q31-qter,  ftp,  lOp.  12ql3-ql5, 
lftp!2-pter,  17q21-qter,  20  (a:  2()p  1 3 )721 7^2 
3ql3-qter.  4pl5,  5,  6pter-ql4  (a:  6p23-pter).  7 
(a:  7q32-qier).  8  (a:  8q22.  8q24).  9pl2-p21.  9q21-qter, 
10q23-qter.  Hql2-q13,  12qT2-q2 1 ,  13q22-qHi% 
14q24-qter,  15q22-qter,  16pI2-pter,  17q,  20  (a:  20pl3. 
20ql3) 

I p34 . 1  -pier.  Iq32-qter.  2ql  1. 2-q  14.1.  5q31-qter,  ftp22-pter, 
8q24,  16pl  1.2-pter,  17q21-qter,  19p,  20,  21,  22qll,2 

lq.  2ql  1 .2-q  14.1.  6p22-pter,  7,  8q24,  9q.  10p,  14q. 
Iftpll.2-pter,  17q21-qter,  ^7^22qlFT 

Iq32-qter.  2ql  1 .2-q  14.1.  5q31-qter,  6p22-pter,  8q24, 
9q21-q22,  lOp,  13q22-qter,  1 4q.  15q22-qtcr. 

1  ftp  1 1 .2-pterT  1 7q2 1  -qter,  20,  21,  22qll.2 
lq32-qter.  2ql  1  -2-q  14.1.  5q31-qter,  6p22-pter,  8q24, 
9q2I-q22,  lOp,  13q22-qter,  14q.  16pl  1.2-pter, 
17q21-qter7~20,  ~21,  22qll.2 
Iq32-qter.  2ql  1 .2-q  14.1.  5q31-qter,  6p22-pter,  8q24, 
9q21-q22,  lOp,  13q22-qter7l4q.  15q22-qter, 

Iftpl  1.2-pter,  17q21-qter,  19p,  20,  21 ,  22qll.2 
Iq32-qter.  2ql  1. 2-q  14.1.  5q31-qter,  ftp22-pter.  8q24, 
9q21-q22,  lOp,  13q22-qter,  14q.  15q22-qter, 

Iftpl  1.2-pter,  17q21-qter,  20,  21,  22q]  1.2 


3p,  4«|,  ftq.  8pter-q21.2,  9p21-pter,  lleen-pl5,  llql4-qter, 
12q21-q23.  Iftq.  18,  19q,  X 

3p.  7cen-q32.  8pter-q21,  8q23.  9q34.  1  lcen-pl 5,  llql4, 

1  lq23-qter.  12q21-qter.  I3een-q22,  18q21-qter,  19q,  X 
3p,  7cen-q32.  8pter-q21,  8q23.  9q34.  1  lcen-pl 5,  llq!4, 

1  lq23-qter.  I2q21-qter.  I3cen-q22.  18q21-qter,  19q,  X 
Icen-q4 1 .  2.  3p.  Ilq23-qter,  13cen-q22.  18q,  X 


3p,  4plft.  4q32-qter,  6ql2-q16,  7p21-q31.  8pter-q21,  8q23. 
9p.  9q32-q34.  IOcen-q25.  Ilcen-pl5,  llql4,  llq23-qter, 
!2p.  1 2q  1 5-qter.  I3cen-q21.  19q,  X 
]pI3-p22.  3p,  4plft.  4q32-qter,  8pter-q21,  8q23,  9p, 
11cen-pl5,  llq)4,  llq23-qter,  12p.  12ql5-qter. 

1 3cen-q2 1 .  Iftq,  18q21-qter,  19q,  X 
3p,  4pl  ft.  4q32-qter,  7q2I-q3L  8pter-q21,  8q23.  9p. 
Ilcen-pl5.  Ilql4,  llq23-qter,  12p.  I2ql5-qter. 
13ccn-q21.  Iftq,  18,  19q7x 
3p  4p  1  ft.  4q32-qter,  7q21-q31.  8pter-q21,  8q23.  9p. 
9q32-q34.  1  lcen-pl  5^11  q  14, Hq2Tqter,  12p, 

1 2q  1 5-qter.  1 3cen-q2 1 .  16q,  18,  19q,  X 
3p,  4plft.  4q32-qter,  7q21-q31.  8pter-q21,  8q23.  9p. 
Ilcen-pl5.  Ilql4.  Ilq23-qter,  1 2p.  12ql 5-qter. 

]  3 cen  q2 1 .  U^TizTX 

3p,  4plft.  4q32-qter»  7q21-q31.  8pter-q21,  8q23.  9p. 

1  lcen-pl 5.  1  lq!4,  llq23-qter,  I2p.  1 2ql 5-qter. 

1 3cen-q2 1 .  iftqTTgpTO^  19q,  X 


a.  Amplil’icaiion 

Bold  and  single  underline,  markers  of  SV40  immortalization 


from  this  process.  The  identification  of  radiation- 
induced  chromosomal  changes  in  the  tumorigenic 
cell  lines  by  use  of  advanced  molecular  cytogenetic 
techniques  such  as  SKY  and  CGH  strongly  sup¬ 
ports  the  usefulness  of  the  established  epithelial 
cell  lines  as  in  vitro  model  for  studying  important 
steps  during  radiation-induced  malignant  transfor¬ 
mation  in  human  thyroid  cells. 
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SUMMARY  High-resolution  physical  maps  are  indispensable  for  directed  sequencing 
projects  or  the  finishing  stages  of  shotgun  sequencing  projects.  These  maps  are  also  critical 
for  the  positional  cloning  of  disease  genes  and  genetic  elements  that  regulate  gene  ex¬ 
pression.  Typically,  physical  maps  are  based  on  ordered  sets  of  large  insert  DNA  clones  from 
cosmid,  P1/PAC/BAC,  or  yeast  artificial  chromosome  (YAC)  libraries.  Recent  technical  devel¬ 
opments  provide  detailed  information  about  overlaps  or  gaps  between  clones  and  pre¬ 
cisely  locate  the  position  of  sequence  tagged  sites  or  expressed  sequences,  and  thus  sup¬ 
port  efforts  to  determine  the  complete  sequence  of  the  human  genome  and  model 
organisms.  Assembly  of  physical  maps  is  greatly  facilitated  by  hybridization  of  non-isotopi- 
cally  labeled  DNA  probes  onto  DNA  molecules  that  were  released  from  interphase  cell  nu¬ 
clei  or  recombinant  DNA  clones,  stretched  to  some  extent  and  then  immobilized  on  a  solid 
support.  The  bound  DNA,  collectively  called  "DNA  fibers,"  may  consist  of  single  DNA  mole¬ 
cules  in  some  experiments  or  bundles  of  chromatin  fibers  in  others.  Once  released  from  the 
interphase  nuclei,  the  DNA  fibers  become  more  accessible  to  probes  and  detection  re¬ 
agents.  Hybridization  efficiency  is  therefore  increased,  allowing  the  detection  of  DNA  tar¬ 
gets  as  small  as  a  few  hundred  base  pairs.  This  review  summarizes  different  approaches  to 
DNA  fiber  mapping  and  discusses  the  detection  sensitivity  and  mapping  accuracy  as  well  as 
recent  achievements  in  mapping  expressed  sequence  tags  and  DNA  replication  sites. 

(J  Histochem  Cytochem  49:939-948,  2001) 
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High-resolution  physical  maps  are  indispensable 
for  large-scale,  cost-effective  gene  discovery.  The  con¬ 
struction  of  such  maps  of  the  human  genome  and 
model  organisms  therefore  has  been  one  of  the  major 
goals  of  the  human  genome  project  (Collins  and  Galas 
1993).  The  most  common  physical  mapping  strategies 
implement  a  bottom-up  approach  for  organizing  indi¬ 
vidual  inserts  from  large  genomic  libraries  into  a  high- 
resolution  map  of  contiguous  overlapping  fragments 
(contigs).  The  progress  in  cloning  large,  megabasepair 
(Mbp)-sized  genomic  DNA  fragments  in  yeast  artifi¬ 
cial  chromosomes  (YACs)  (Cohen  et  al.  1993;  Olson 
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1993)  has  made  it  possible  to  rapidly  construct  low- 
resolution  framework  maps  based  on  overlapping 
YAC  clones.  Fluorescence  in  situ  hybridization  (FISH) 
has  proved  indispensable  for  identification  of  non-chi- 
meric  YAC  clones  and  for  physical  mapping  of  indi¬ 
vidual  YAC  clones  onto  metaphase  chromosomes 
(Selleri  et  al.  1992).  To  generate  physical  maps  with  a 
resolution  of  1  Mbp  or  better,  the  YAC  clones  are  or¬ 
dered  by  combining  different  complementing  analytic 
techniques:  pulsed-field  gel  electrophoresis  (PFGE) 
(Vetrie  et  al.  1993),  FISH  with  interphase  cell  nuclei 
or  metaphase  spreads  (Trask  et  al.  1989;  Brandriff  et 
al.  1992;  Ashworth  et  al.  1995),  sequence  tagged  site 
(STS)  content  mapping  (Green  and  Olson  1990;  Cof¬ 
fey  et  al.  1992;  Weissenbach  et  al.  1992;  Lu-Kuo  et  al. 
1994;  Locke  et  al.  1996),  and/or  DNA  repeat  finger¬ 
printing  (Bellanne-Chantelot  et  al.  1992;  Zucchi  and 
Schlessinger  1992;  Waterston  and  Sulston  1995;  Bell 
et  al.  1995). 
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High-resolution  maps  providing  ordered  sets  of 
cloned  DNA  fragments  at  the  100-kb  level  of  resolu¬ 
tion  are  assembled  with  smaller,  more  manageable 
DNA  fragments  isolated  from  other  libraries.  Most 
groups  prefer  cloning  of  genomic  DNA  in  vectors  that 
maintain  relatively  large  DNA  fragments  without  re¬ 
arrangements,  are  non-chimeric,  and  allow  easy  DNA 
purification.  In  general,  high-resolution  maps  are 
composed  of  overlapping  cosmids  (Stallings  et  al. 
1992;  Tynan  et  al.  1992;  Nizetic  et  al.  1994),  Pl/PAC 
clones  (Pierce  et  al.  1992;  Ioannou  et  al.  1994),  or 
bacterial  artificial  chromosomes  (BACs)  (Shizuya  et  al. 
1992).  Assembly  of  the  high-resolution  maps  requires 
identification  of  cloned  DNA  sequences  that  contain 
overlapping  regions  of  the  genome  (Cheng  and  Weier 
1997).  To  minimize  the  effort  invested  in  assembling  a 
contiguous  set  of  clones,  clone  overlaps  need  to  be  de¬ 
termined  quickly  and  accurately.  This  has  been  ac¬ 
complished  by  various  forms  of  clone  fingerprinting 
[e.g.,  by  identification  of  common  restriction  fragment 
or  inter-Alu  PCR  patterns  (Branscomb  et  al.  1990; 
Nelson  1991;  Patil  et  al.  1994),  by  hybridization  to 
bound  filter  clone  arrays  (Hoheisel  and  Lehrach  1993; 
Locke  et  al.  1996)  or  oligonucleotide  arrays  (Sapolsky 
and  Lipshutz  1996),  and  by  identification  of  overlap¬ 
ping  sequence  tagged  sites  (Green  et  al.  1991;  Abura- 
tani  et  al.  1996).  The  development  of  radiation  hybrid 
mapping  (Cox  et  al.  1990),  the  availability  of  large 
numbers  of  STS  markers,  and  the  construction  of  ex¬ 
tensive  bacterial  clone  resources  provide  additional 
means  to  accelerate  the  process  of  mapping  a  chromo¬ 
some  and  preparing  clone  contigs  ready  for  sequenc¬ 
ing.  These  techniques,  although  effectively  used  by  the 
genome  community,  are  limited  because  they  do  not 
readily  yield  information  about  contig  orientation,  ex¬ 
tent  of  deletions  or  rearrangements  in  clones,  overlap 
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of  contig  elements,  or  their  chimerism  status,  nor  do 
they  provide  information  about  the  extent  of  gaps  in 
the  maps.  The  precise  localization  of  cloned  DNA 
fragments  within  much  larger  genomic  fragments  and 
knowledge  about  the  extent  of  overlap  between  two 
clones  is  needed  to  assemble  high-resolution  physical 
maps.  Techniques  to  rapidly  identify  minimally  over¬ 
lapping  clones  and  to  determine  the  extent  and  orien¬ 
tation  of  overlap  expedite  the  construction  of  minimal 
tiling  paths  and  facilitate  the  sequence  assembly  pro¬ 
cess.  As  demonstrated  in  this  review,  FISH  can  pro¬ 
vide  this  critical  information.  Furthermore,  FISH- 
based  mapping  is  one  of  the  few  mapping  procedures 
without  requirements  for  different  clones  to  overlap  to 
be  positioned  and  ordered  (Horelli-Kuitunen  et  al. 
1999). 

In  early  applications  of  FISH-based  clone  ordering 
and  assembly  of  physical  maps,  Brandriff  and  col¬ 
leagues  (1991a)  took  advantage  of  the  high  hybridiza¬ 
tion  efficiency  and  enhanced  spatial  resolution  ob¬ 
tained  with  decondensed  human  sperm  DNA.  In  their 
approach,  individual  human  sperm  nuclei  were  fused 
with  hamster  eggs  and  probes  were  hybridized  to  the 
decondensed  DNA  of  human  pronuclei  (Brandriff  et 
al.  1991b).  The  example  in  Figure  IB  demonstrates 
such  decondensation  in  form  of  a  cloud  or  network  of 
green  signals  produced  by  a  biotin-labeled  probe  di¬ 
rected  towards  the  repeated  satellite  III  DNA  in  the 
heterochromatic  region  of  chromosome  9qll.  The 
probe  visualized  with  avidin-FITC  fluoresced  in  green 
(Figure  IB),  and  the  DNA  was  counterstained  with 
DAPI  emitting  blue  fluorescence  (Figure  1A).  The  hu¬ 
man  sperm-hamster  egg  fusion  technique  was  used 
successfully  to  order  cosmid  clones  in  preparing  high- 
resolution  physical  maps  for  human  chromosome  19 
(Brandriff  et  al.  1992;  Ashworth  et  al.  199 5). 


Figure  1  DNA  fiber  mapping.  (A,B)  Hybridization  to  decondensed  sperm  nuclei  DNA  allows  probes  to  be  ordered  with  high  resolution.  (B) 
Results  of  hybridizing  a  biotinylated  probe  specific  for  the  chromosome  9-specific  satellite  III  DNA  and  detection  of  bound  probe  molecules 
with  avidin-FITC.  (A)  DAPI-counterstained  DNA.  (C)  DNA  released  from  cell  nuclei  can  be  hybridized  very  efficiently  with  non-isotopically  la¬ 
beled  probes.  A  biotinylated  DNA  probe  specific  for  a-satel I ite  repeats  on  chromosomes  13  and  21  generated  hybridization  signals  that  ap¬ 
pear  yellow  due  to  overlap  between  the  green  probe  fluorescence  and  the  red  fluorescence  of  propidium  iodide,  which  was  used  to  coun¬ 
terstain  the  DNA.  (D-G)  High-resolution  mapping  of  a-satellite  DNA  domains  using  surface  spreads  of  DNA  fibers.  Unfixed  human  white 
blood  cells  were  spread  on  glass  slides,  fixed  with  acetic  acid:methanol,  and  hybridized  with  a  biotinylated  probe  for  chromosome  13/21- 
specific  alphoid  DNA.  Green  signal  (D,F)  delineating  the  alphoid  tandem  repeats  appear  as  beads-on-a-string  in  the  more  extended  chroma¬ 
tin  parts.  Images  of  the  propidium  iodide-counterstained  DNA  (E,G)  demonstrate  hybridization  tracks  located  in  DNA  fibers  of  very  differ¬ 
ent  size.  (H)  Lambda  DNA  molecules  bound  to  an  APS-derivatized  glass  slide  were  straightened  by  the  receding  meniscus  during  drying  and 
hybridized  with  biotin-  or  digoxigenin-labeled  restriction  fragments  which  were  detected  with  avidin-FITC  (green)  and  rhodamine-conju- 
gated  antibodies  against  digoxigenin  (red).  (I)  Circular  DNA  molecules  were  excised  from  a  PFGE  gel  and  purified  by  agarose  digestion.  A 
small  plasmid  probe  of  about  2  kb  was  mapped  onto  a  single,  closed  circular  BAC  DNA  molecule.  An  FITC-labeled  probe  mix  (green)  delin¬ 
eates  the  BAC  vector  A  small  probe  (-1200  bp)  binds  within  the  vector  part  near  the  T7  promoter  (red).  The  BAC  molecule  was  counter- 
stained  by  hybridization  of  a  probe  detected  with  avidin  conjugated  to  AMCA  (blue).  (J,K)  The  extent  of  overlap  between  clones  can  be  de¬ 
termined  rapidly  by  hybridizing  one  clone  onto  DNA  fibers  prepared  from  the  other  clone.  Here,  a  BAC  clone  (red)  was  mapped  onto  a 
partially  overlapping  circular  PAC  molecule  (blue).  The  vector  part  of  about  16  kb  was  counterstained  with  green/red  probe  combination. 
The  results  are  summarized  schematically  in  K.  (L)  Hybridization  of  labeled  cDNA  (red  signals  on  top  of  blue  signals)  onto  DNA  fibers  pre¬ 
pared  from  genomic  clones  such  as  BACs  reveals  the  position  and  genomic  organization  of  expressed  sequences.  The  BAC  molecule  was 
counterstained  by  hybridization  of  a  probe  detected  with  avidin  conjugated  to  AMCA  (blue).  The  vector  part  of  about  8  kb  was  counter¬ 
stained  with  green/red  probe  combination.  Bars  =  10  |j.m. 
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The  procedure,  however,  had  two  obvious  prob¬ 
lems.  Fusion  of  human  sperm  with  hamster  eggs  and 
fixation  of  pronuclei  is  a  time-consuming,  laborious 
process  and  might  not  scale  well  enough  to  meet 
the  high-throughput  requirements  of  most  genome 
projects.  The  second  shortfall  of  the  procedure  was  a 
complete  lack  of  control  over  the  extent  of  DNA 
decondensation  and  orientation  of  pronuclei.  This 
spurred  efforts  in  the  1990s  to  manipulate  chromatin 
or  purified  DNA  molecules  that  could  serve  as  a  tem¬ 
plate  for  high-resolution  physical  mapping  of  DNA 
probes.  The  optimal  procedure  would  be  inexpensive, 
rapid,  reproducible,  and  deliver  mapping  data  limited 
only  by  the  resolution  of  the  light  microscope.  A  de¬ 
cade  later,  we  find  ourselves  equipped  with  an  arsenal 
of  complementary  FISH-based  mapping  procedures 
that  cover  a  very  broad  range  of  mapping  intervals. 
Furthermore,  the  simultaneous  development  of  more 
sensitive  fluorescence  detection  reagents  has  pushed 
the  limits  of  detection  down  to  a  few  hundred  base- 
pairs  (bp). 

DNA  Fiber  Mapping 

The  expression  “DNA  fiber  mapping”  has  become  a 
collective  name  for  quite  different  mapping  tech¬ 
niques.  As  indicated  in  Figure  2,  the  diameter  of 
DNA  fibers  increases  as  DNA  molecules  with  a  di¬ 
ameter  of  2  nm  are  packed  into  chromatin  ranging 
from  10  nm  for  histone-packed  DNA  molecules  and 
30  nm  for  chromatin  fibers  all  the  way  to  chromatids 
of  700  nm  diameter.  Chromatin  can  be  released  from 
interphase  cell  nuclei  by  various  chemical  or  mechani¬ 
cal  methods,  and  investigators  tried  to  coin  names  that 
reflect  their  individual  approach.  Isolation  of  DNA 
from  cell  nuclei,  extension,  and  preparation  of  chro¬ 
matin  or  DNA  fibers  with  diameters  ranging  in  size 
from  a  few  to  several  hundred  nm  (Figure  2)  improves 
the  accessibility  of  the  DNA  targets  for  both  probes 
and  detection  reagents.  Accordingly,  the  hybridization 
efficiencies  increase,  and  DNA  targets  of  less  than  1 
kb  can  be  detected  routinely  using  procedures  nor¬ 
mally  applied  in  metaphase  and  interphase  cell  FISH. 

In  1992,  Heng  et  al.  described  the  use  of  chemicals 
to  release  chromatin  from  interphase  cell  nuclei.  The 
results  look  somewhat  similar  to  the  propidium  io¬ 
dide-stained  free  chromatin  shown  in  Figure  1C.  Heng 
et  al.  (1992)  used  the  drug  N-[4-(9-acridinylamino)- 
3-methoxyphenyl]-methanesulfonamide  (ra-AMSA)  and 
an  alkaline  lysis  procedure  to  release  the  DNA  from 
interphase  cells  and  demonstrated  mapping  of  DNA 
sequences  in  distances  ranging  from  21  kb  to  350  kb 
onto  the  extended  chromatin.  The  authors  speculated 
that,  with  multicolor  FISH,  free  chromatin  mapping 
would  readily  resolve  gene  sequences  separated  by  as 
little  as  10  kb  (Heng  et  al.  1992).  Soon  thereafter,  Sen- 
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ger  et  al.  (1994)  and  Fidlerova  et  al.  (1994)  demon¬ 
strated  the  use  of  a  sodium  hydroxide/ethanol  mixture 
or  70%  formamide  to  prepare  chromatin  fibers  which 
they  called  “free  DNA.”  Hybridization  signals  from 
cosmid  probes  appeared  as  extended  lines  and  over¬ 
laps  between  clones  were  determined  using  the  lengths 
of  the  probe  signals  and  their  known  kilobase  size  as 
an  internal  standard.  Typical  experiments  calculated 
averages  from  30  or  more  signals  (Senger  et  al.  1994). 

Wiegant  et  al.  (1992)  and  Lawrence  et  al.  (1992) 
applied  a  different  chemical  method  to  release  DNA 
from  cells  producing  so-called  “DNA  halos.”  The  pro¬ 
cedure  was  based  on  detergent  and  high  salt  ex¬ 
traction  followed  by  intercalating  dye/UV  treatment 
(Cook  et  al.  1976;  Vogelstein  et  al.  1980).  Using  mul¬ 
ticolor  hybridizations,  these  investigators  were  able  to 
determine  relative  map  positions  and  to  detect  10-kb 
overlap  between  individual  cosmid  clones.  Wiegant  et 
al.  (1992)  noted  cosmid-specific  signals  that  appeared 
as  linear  beaded  signals  of  about  10  |xm  (i.e.,  about 
3.5-4  kb/fxm),  and  the  authors  concluded  that  their 
DNA  halo  preparations  produced  “essentially  linear¬ 
ized”  DNA.  The  map  configuration  of  a-satellite 
DNA  arrays  or  single-copy  DNA  could  typically  be 
derived  from  analysis  of  only  5-10  cells.  These  find¬ 
ings  were  confirmed  independently  by  the  work  of 
Lawrence  et  al.  (1992).  But  as  Gerdes  et  al.  pointed 
out  in  their  publication  in  1994,  “many  procedures 
can  be  used  to  decondense  DNA  for  in  situ  hybridiza¬ 
tion,  however  many  of  these  result  in  highly  variable 
and  apparently  random  threads  of  DNA  smeared  er¬ 
ratically  across  the  slide,  with  little  if  any  retention  of 
native  morphology”  (Gerdes  et  al.  1994).  As  the  im¬ 
ages  in  Figures  1D-1G  illustrate,  DNA  halo  prepara¬ 
tions  often  show  residual  nuclei  containing  compact 
hybridization  domains,  while  the  direction,  extent, 
and  straightness  of  hybridization  tracks  on  the  ex¬ 
tended  chromatin  appear  to  be  highly  variable. 

If  DNA  molecules  could  be  stretched  uniformly  in 
one  direction,  they  might  provide  linear  templates  for 
visual  FISH  mapping.  Parra  and  Windle  (1993)  de¬ 
scribed  a  procedure  they  termed  “direct  visual  hybrid¬ 
ization”  (DIRVISH),  in  which  100-5000  cells  in  phos¬ 
phate  buffer  were  placed  on  a  one  end  of  a  slide  and 
air-dried.  The  slide  was  then  immersed  in  a  lysis  solu¬ 
tion  containing  detergent,  removed  after  5  min,  and 
tilted  to  allow  the  drop  of  DNA  to  run  down  the  slide. 
The  procedure  produced  some  linear  arrangements  of 
DNA  but,  on  the  basis  of  the  extent  of  hybridization 
signals,  the  authors  estimated  that  the  DNA  was  ex¬ 
tended  to  a  length  ranging  from  40%  to  200%  of  the 
length  of  relaxed  duplex  DNA,  i.e.,  3.4  nm/bp  (Parra 
and  Windle  1993). 

In  the  following  year,  Heiskanen  et  al.  (1994)  pub¬ 
lished  a  protocol  to  stretch  genomic  DNA  embedded 
in  agarose.  The  DNA  was  melted  on  a  microscope 
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Figure  2  Packaging  of  DNA  into 
chromatin  and  chromosomes.  Several 
levels  of  packing  of  DNA  can  be  dis¬ 
criminated.  Two-nanometer-wide  DNA 
molecules  are  folded  into  1 1-nm  chro¬ 
matin  segments  which,  when  tightly 
packed  with  nudeosomes,  appear  as 
30-nm  chromatin  fibers.  Further  pack¬ 
ing  leads  to  chromosomal  sections 
with  a  diameter  of  about  300  nm  and 
the  700-nm-thick  chromatids  of 
metaphase  chromosomes. 


slide  by  heating  the  slide,  then  distributed  and 
stretched  manually  using  a  coverslip.  Major  draw¬ 
backs  of  this  manual  procedure  were  heterogeneous 
stretching  and  a  rather  low  abundance  of  informative 
signals  on  the  slides. 

Procedures  published  by  1995  allowed  FISH  for 
most  preparations  of  decondensed  nuclear  or  isolated 
cloned  DNA  and  visualization  of  probe  overlap  to 
provide  some  information  about  the  existence  and  size 
of  gaps  between  clones  (Parra  and  Windle  1993;  Flor- 
ijn  et  al.  1995;  Heiskanen  et  al.  1995,1996).  However, 
none  of  these  techniques  provided  sufficiently  accurate 
information  about  the  extent  of  clone  overlap  or  the 
separation  between  elements  in  the  map  because  the 
chromatin  onto  which  clones  were  mapped  was  con¬ 
densed  to  various  degrees  from  site  to  site  in  these 
preparations  (Rosenberg  et  al.  1995). 

We  demonstrated  that  cloned  DNA  fragments  can 
readily  be  mapped  by  FISH  onto  DNA  molecules 
straightened  by  the  hydrodynamic  action  of  a  receding 


meniscus  and,  referring  to  its  quantitative  nature,  we 
termed  our  technique  “Quantitative  DNA  Fiber  Map¬ 
ping  (QDFM)”  (Weier  et  al.  1995).  In  QDFM,  a  solu¬ 
tion  of  purified  DNA  molecules  is  placed  on  a  flat  sur¬ 
face  prepared  so  that  the  DNA  molecules  slowly 
attach  at  one  or  both  ends.  The  DNA  solution  is  then 
spread  over  a  larger  area  by  placing  a  coverslip  on  top. 
DNA  molecules  are  allowed  to  bind  to  the  surface. 
During  drying,  the  molecules  are  straightened  and  uni¬ 
formly  stretched  by  the  hydrodynamic  action  of  the 
receding  meniscus.  Molecules  prepared  in  this  manner 
are  stretched  with  remarkable  homogeneity.  We  esti¬ 
mated  that  a  properly  stretched  molecule  should  ex¬ 
tend  about  ^2.3  kb/jjim,  i.e.,  approximately  30%  over 
the  length  predicted  for  a  double-stranded  DNA  mole¬ 
cule  of  the  same  size  (Bensimon  et  al.  1994;  Weier  et 
al.  1995)  (Figure  1H).  We  also  showed  that  QDFM 
can  be  applied  to  linear  or  circular  DNA  molecules 
ranging  in  size  from  a  few  kb  to  more  than  1  Mbp, 
which  enabled  us  to  determine  clone  overlaps  or  to 
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map  small  probes  with  near-kilobase  resolution  onto 
entire  yeast  chromosomes  and  mega- Y AC  clones  from 
the  CEPH/Genethon  library  (Wang  et  al.  1996)  (Fig¬ 
ures  1J— IK). 

The  subsequent  studies  of  Hu  et  al.  (1996)  and 
Yokota  et  al.  (1997)  presented  rapid  protocols  to  pre¬ 
pare  DNA  fibers  on  3-aminopropyl-triethoxysilane- 
coated  glass  or  mica  slides.  In  their  studies  imaging  X 
DNA  molecules  by  atomic  force  microscopy,  Hu  et  al. 
(1996)  let  coverslipped  DNA  molecules  bind  to  the 
surface  for  only  a  few  minutes.  When  the  coverslip 
was  removed  for  imaging,  the  molecules  were  straight¬ 
ened  by  the  moving  meniscus.  Yokota  et  al.  (1997)  de¬ 
veloped  a  protocol  in  which  the  meniscus  motion  was 
controlled  mechanically,  which  provided  advantages 
in  speed  and  uniformity  of  the  straightened  molecules. 

Overcoming  Challenges  in  Mapping  Tandemly 
Repeated  DNA  Sequences 

Most  mammalian  genomes  contain  very  large  blocks 
of  heterochromatin.  Certain  regions  such  as  the 
(near)-centromeric  heterochromatin  seem  to  be  in¬ 
volved  in  karyokinesis  and  chromosome  association 
during  meiosis  (Dozortsev  et  al.  2000),  but  the  overall 
role  of  these  clusters  of  “junk  DNA”  remains  largely 
unknown  (Epplen  et  al.  1998).  Physical  map  assembly 
in  heterochromatic  regions  by  more  conventional 
methods,  i.e.,  clone  fingerprinting  or  STS  content 
mapping,  has  proved  less  efficient,  and  initial  releases 
of  maps  of  the  human,  Drosophila,  and  other  complex 
genomes  might  suffer  from  the  presence  of  major 

gaps- 

Several  groups  applied  DNA  fiber  mapping  tech¬ 
niques  to  investigate  the  interface  between  euchro- 
matic  and  heterochromatic  regions  as  well  as  the  orga¬ 
nization  of  tandem  DNA  repeats  in  humans  (Wiegant 
et  al.  1992;  Haaf  and  Ward  1994b;  Shiels  et  al.  1997), 
animals  (Sjoberg  et  al.  1997;  Li  et  al.  2000),  and 
plants  (Fransz  et  al.  1996;  Jackson  et  al  1998,1 999; 
Zhong  et  al.  1998).  The  methods  applied  by  these  re¬ 
searchers  ranged  from  DNA  halo  preparations  and  ex¬ 
tended  chromatin  to  highly  purified  DNA  molecules. 
Published  results  showed  very  specific  beads-on-a-string 
patterns  and  presented  convincing  cases  that  DNA  fi¬ 
ber  mapping  is  a  powerful  tool  to  decipher  the  organi¬ 
zation  and  extent  of  centromeric  and  telomeric  DNA 
repeats. 

Mapping  of  Single-copy  DNA 

Physical  Map  Assembly 

The  assembly  of  high-resolution  physical  maps  has 
been  a  major  application  of  DNA  fiber  mapping.  Suc¬ 
cessful  efforts  to  construct  local  maps  for  delineation 
of  disease  genes  were  mentioned  above,  but  the  suc¬ 


cess  of  additional  investigators  deserves  proper  recog¬ 
nition. 

Klockars  et  al.  (1996)  used  fiber-FISH  to  construct 
a  physical  map  of  the  CLN5  region  on  13q22.  Duell  et 
al.  (1997,1998)  constructed  high-resolution  physical 
maps  based  on  QDFM  for  the  immunoglobulin-(3  re¬ 
gion  on  chromosome  22  and  a  large  region  including 
the  5'  end  of  the  human  apolipoprotein  B  gene  on 
chromosome  2.  Theuns  et  al.  (1999)  determined  the 
genomic  organization  of  the  human  presenilin  1  gene, 
which  was  localized  on  chromosome  14q24.3,  by  fi¬ 
ber-FISH.  Although  applying  different  techniques  to 
prepare  the  DNA  fibers,  these  recent  applications  of 
DNA  fiber  mapping  were  able  to  support  ongoing 
mapping  efforts  by  providing  physical  maps  with 
near-kilobase  resolution. 

Direct  Visualization  of  Gene  Amplifications 
and  Deletions 

Linearly  extended  chromatin  or  DNA  molecules  are 
ideal  substrates  to  study  gene  amplification  or  dele¬ 
tions.  Parra  and  Windle  (1993)  applied  the  DIRVISH 
technique  to  resolve  low  level  amplifications  of  the  di¬ 
hydrofolate  reductase  (DHFR)  gene  in  a  hamster  cell 
line.  All  DNA  fiber  mapping  procedures  are  single 
DNA  molecule  investigations,  which  permit  study  of 
cell-to-cell  variations.  Parra  and  Windle  (1993)  pre¬ 
sented  evidence  that  DIRVISH  can  resolve  the  diver¬ 
sity  in  structure  and  complexity  found  in  the  amplified 
DHFR  array,  which  would  normally  escape  detection 
and  proper  characterization. 

Michalet  et  al.  (1997)  extended  human  genomic 
DNA  on  silanized  glass  surfaces  using  a  receding  me¬ 
niscus  and  found  that  a  precise  measurement  of  hy¬ 
bridized  DNA  probes  can  be  achieved  without  nor¬ 
malization.  Stretching  yeast  genomic  DNA  amounts 
representing  hundreds  of  eukaryotic  genomes,  the  au¬ 
thors  performed  the  high-resolution  mapping  of 
cosmid  contigs  on  a  yeast  artificial  chromosome 
(YAC).  Furthermore,  stretching  of  human  genomic 
DNA  allowed  measurement  of  gaps  between  contigs 
and  determination  of  microdeletions  in  the  tuberous 
sclerosis  2  gene  of  patients’  DNA. 

Mapping  of  Expressed  Sequences 

Applications  of  DNA  fiber  mapping  extend  beyond 
map  assembly  and  can  provide  valuable  information 
for  clone  validation,  definition  of  a  minimal  tiling 
path,  and  quality  control  in  the  sequence  assembly 
process.  Even  more  exciting,  its  high  hybridization  ef¬ 
ficiency  makes  DNA  fiber  mapping  the  method  of 
choice  for  visual  mapping  of  expressed  sequences. 

Several  approaches  exist  to  map  expressed  se¬ 
quences  or  to  study  the  organization  of  larger  genes. 
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Gene  fragments  or  entire  cDNAs  cloned  in  plasmid 
vectors  can  be  amplified  by  in  vitro  DNA  amplifica¬ 
tion  using  the  polymerase  chain  reaction  (PCR),  la¬ 
beled  with  reporter  molecules,  and  used  as  hybridiza¬ 
tion  probes  (Wang  et  al.  1996).  Figure  1  shows  two 
examples  of  this  application.  A  cloned  cDNA  frag¬ 
ment  containing  exon  2  of  the  human  protein  4.1  gene 
(Peters  et  al.  1998)  (Figure  II)  was  positioned  accu¬ 
rately  on  a  closed  circular,  co-linear  BAC  DNA  mole¬ 
cule.  Counterstaining  the  BAC  vector  part  with  a  com¬ 
bination  of  FITC-  and  digoxigenin-labeled  probes  that 
were  detected  in  green  and  red,  respectively,  the  posi¬ 
tion  of  gene  4.1  exon  2  within  the  BAC  insert  (the  red 
signal  on  top  of  blue  counterstain)  can  be  determined 
with  near-kilobase  resolution. 

The  organization  of  genes  in  transcribed  sequences 
interrupted  by  intronic  sequences  can  be  easily  dem¬ 
onstrated  by  DNA  fiber-FISH.  We  used  QDFM  in  an 
attempt  to  resolve  the  structure  of  a  gene  from  chro¬ 
mosome  20  frequently  found  amplified  in  human  tu¬ 
mors  (Tanner  et  al.  1994).  As  Figure  1L  shows,  a 
cDNA  clone  of  about  5  kb  maps  to  three  distinct  posi¬ 
tions  on  a  co-linear  BAC  clone.  This  result  should  be 
considered  a  rough  estimate  because  exons  smaller 
than  500  bp  may  not  be  detectable  with  this  simple 
method  (Weier  et  al.  1995).  On  the  other  hand, 
QDFM  mapping  turned  out  to  be  a  highly  reproduc¬ 
ible  method,  so  that  was  sufficient  to  analyze  no  more 
than  five  to  eight  BAC  molecules  (Weier  2001). 

Technical  developments  in  recent  years  have 
opened  the  doors  to  mapping  of  even  smaller  exons. 
In  1996,  Florijn  et  al.  demonstrated  the  use  of  DNA 
halo  preps  for  mapping  of  exons  ranging  in  size  from 
202  to  778  bp.  Co-hybridization  of  co-linear  cosmid 
clones  enabled  this  group  to  reproducibly  locate  exon 
fragments  of  about  200  bp  on  extended  genomic  DNA 
in  the  context  of  the  cognate  cosmid  signal.  Detection 
efficiencies  of  70-90%  were  found  with  probes  larger 
than  400  bp,  but  the  detection  efficiency  decreased  to 
about  30%  when  fragments  of  about  200-250  bp 
were  mapped. 

More  recently,  Aaltonen  and  co-workers  (1997; 
Horelli-Kuitunen  et  al.  1999)  used  a  tyramide-based 
method  (Raap  et  al.  1995)  for  detection  of  small  hy¬ 
bridization  targets.  By  applying  this  sensitive  detection 
scheme  to  create  a  visual  transcript  map  of  six  genes 
on  chromosome  21q22.3,  Horelli-Kuitunen  et  al. 
(1999)  were  able  to  detect  targets  ranging  from  316  to 
956  bp  on  25-44%  of  the  fibers.  Independent  experi¬ 
ments  confirmed  that  an  average  of  89%  of  the  EST 
hybridization  signals  were  specific  signals.  These  re¬ 
sults  suggested  that  a  mere  increase  in  detection  sensi¬ 
tivity  may  not  be  sufficient  to  significantly  increase  the 
success  rates  of  DNA  fiber  mapping  experiments,  i.e., 
the  overall  efficiency  of  hybridization  and  detection. 

In  the  experiments  conducted  by  Horelli-Kuitunen 


et  al.  (1999),  the  DNA  fibers  were  prepared  from  aga¬ 
rose-embedded  human  genomic  lymphocyte  chroma¬ 
tin  according  to  the  protocol  of  Heiskanen  et  al. 
(1994,1995),  which  may  have  limited  the  access  of 
probes  and  detection  reagents  to  the  target  DNA.  As 
van  der  Rijke  et  al.  (2000)  proposed,  the  loss  of  target 
DNA  and  limited  accessibility  due  to  in  situ  renatur- 
ation  and  attachment  might  have  a  negative  impact  on 
the  detection  efficiency.  In  the  author’s  opinion,  DNA 
fiber-FISH  sensitivity  is  determined  by  hybridization 
efficiency,  rather  than  the  ability  to  generate  sufficient 
signal  from  small  probes  (van  der  Rijke  et  al.  2000). 

High-resolution  Studies  of  DNA  Replication 

Despite  intense  efforts,  the  orderly  activation  of  repli¬ 
cation  sites  in  genomes  of  higher  organisms  remains 
largely  unexplained.  The  main  reason  for  this  may  be 
the  complexity  of  a  process  orchestrating  the  partly 
parallel  activation  of  an  estimated  104  to  106  replica¬ 
tion  sites. 

DNA  fiber  mapping,  with  its  high  resolution  and 
sensitivity,  might  provide  important  mapping  infor¬ 
mation  about  the  location  and  spacing  of  replication 
sites.  Several  groups  used  fiber  mapping  techniques  to 
demonstrate  replication  forks  in  yeast  (Rosenberg  et 
al.  1995),  in  Xenopus  laevis  (Herrick  et  al.  2000),  and 
in  humans  (Gerdes  et  al.  1994;  Haaf  1996).  The  ap¬ 
proach  chosen  by  Rosenberg  et  al.  (1995)  was  based 
on  the  previously  described  DIRVISH  and  DNA  halo 
techniques  (Wiegant  et  al.  1992;  Parra  and  Windle 
1993).  Hybridizing  differentially  labeled  cosmid 
probes  to  yeast  and  YAC  DNA  fibers,  the  authors 
could  provide  convincing  visual  evidence  of  DNA  rep¬ 
lication  “forks”  and  “eyes”  in  YACs  containing  hu¬ 
man  DNA  (Rosenberg  et  al.  1995).  A  rather  limited 
number  of  experiments  enabled  the  authors  to  identify 
two  origins  of  replication  in  a  400-kb  region  contain¬ 
ing  the  human  dystrophin  gene,  to  present  evidence 
for  bi-directional  replication  and  to  demonstrate  that 
human  DNA  cloned  in  yeast  is  capable  of  initiating  its 
own  replication. 

Studies  described  by  Gerdes  et  al.  (1994),  Haaf 
(1996),  and  Herrick  et  al.  (2000)  used  incorporation 
of  non-isotopically  labeled  nucleotides  as  markers  for 
replicated  DNA.  The  results  of  Gerdes  et  al.  (1994) 
showed  pulse-incorporated  bromodeoxyuridine  (BrdU) 
representing  nascent  replicating  DNA  localized  with 
the  base  of  the  chromatin  loops  in  halo  preparations. 
At  increasing  chase  times,  the  replicated  DNA  was 
consistently  found  farther  out  on  the  extended  region 
of  the  halo  (Gerdes  et  al.  1994).  Haaf  (1996)  investi¬ 
gated  the  timing  of  replication  of  a-satellite  DNA  ar¬ 
rays  by  incorporation  of  BrdU  into  newly  synthesized 
DNA,  followed  by  hybridization  of  biotinylated  chro¬ 
mosome-specific  alphoid  probes  to  DNA  released 
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from  fibroblast  or  lymphoblast  cell  nuclei.  Using  the 
DNA  halo  preparations,  he  was  able  to  demonstrate 
many  independent  replication  sites  within  the  large 
clusters  of  tandemly  repeated  alphoid  DNA  and  par¬ 
tially  hemizygous  hybridization  patterns  in  support  of 
the  hypothesis  that  replication  of  a-satellite  DNA  on 
homologous  chromosomes  is  highly  asynchronous. 

The  recent  work  of  Herrick  et  al.  (2000)  extended 
previous  applications  of  DNA  fibers  produced  by  the 
hydrodynamic  action  of  a  receding  meniscus  (Bensi- 
mon  et  al.  1994;  Weier  et  al.  1995).  Xenopus  sperm 
nuclei  were  labeled  by  incorporation  of  biotin-dUTP 
in  newly  synthesized  DNA  followed  by  pulse-labeling 
with  digoxigenin-dUTP  at  fixed  time  intervals.  Be¬ 
cause  triphosphate  deoxynucleosides  do  not  efficiently 
traverse  the  intact  cell  membrane,  all  experiments  had 
to  be  carried  out  with  de-membranated  sperm  nuclei 
and  activated  egg  extracts  (Hyrien  and  Mechali  1992). 
Using  red  and  green  fluorescent  antibodies  directed 
against  biotin  and  digoxigenin,  respectively,  the  au¬ 
thors  were  able  to  discriminate  early  from  late  repli¬ 
cating  sequences  and  to  obtain  estimates  of  the  forma¬ 
tion  of  new  replication  forks. 

Conclusions 

Recent  research  has  lead  to  major  improvements  in 
hybridization-based  physical  mapping  procedures.  De¬ 
velopments  in  DNA  fiber  mapping  have  reached  the 
point  of  immediate  practical  utility:  DNA  probes  as 
small  as  500  bp  can  be  mapped  routinely  onto  immo¬ 
bilized  templates  composed  of  linear  or  circular  DNA 
molecules  that  range  in  size  from  less  than  10  kb  to 
more  than  one  Mbp.  Uniform  stretching  facilitates  the 
conversion  of  measured  physical  distances  into  ge¬ 
nomic  distances.  The  highly  reproducible  stretching 
procedures  require  analysis  of  only  a  few  DNA  mole¬ 
cules  for  accurate  determination  of  map  positions  by 
multicolor  fluorescence  microscopy  and  digital  image 
analysis.  A  mapping  accuracy  in  the  kb  range  coupled 
to  efficient  signal  amplification  procedures  to  visualize 
signals  from  small  targets  allows  rapid  assembly  of 
high-resolution  physical  maps  for  large-scale  sequenc¬ 
ing  and  map  closure  as  well  as  high-resolution  maps 
of  expressed  sequences. 

DNA  fiber  mapping  technologies  will  enhance  per¬ 
formance  of  virtually  all  mapping  and  sequencing 
projects,  including  ongoing  and  future  sequencing  of 
model  organisms  and  bacterial  genomes.  Implementa¬ 
tion  of  these  technologies  will  expedite  sequencing  by 
increasing  parallelism,  will  lower  the  overall  cost  by 
reducing  template  redundancies,  and  will  expedite 
map  closure.  Furthermore,  the  techniques  will  benefit 
the  positional  cloning  of  disease  genes  and  character¬ 
ization  of  genomic  elements  controlling  their  expres¬ 
sion. 
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